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Abstract. We examine the fluctuations of the empirical density measure for the 
colour version of the symmetric nearest neighbour zero range particle systems in 
dimension one. We show that the weak limit of these fluctuations is the solution of 
a system of coupled generalized Ornstein-Uhlenbeck processes. We also discuss how 
this result may be used to prove a central limit theorem for the tagged particle on 
the level of finite dimensional distributions, and identify the limiting variance. This 
is the central limit theorem associated to propagation of chaos for this interacting 
particle system. 

1. Introduction 

The zero range particle system describes a class of microscopic models. It was 
first introduced by Spitzer in 1970 |Spi| as an example of an interacting particle 
system, and has since been studied at length. One of its greatest advantages is its 
mathematical tractability. It has found widespread application in the modelling of 
nonequilibrium phenomena. Examples of these include models of sandpile dynamics 
and other flow mechanisms, as well as the repton model of gel electrophoresis |Evaj . 

In the zero range system we consider particles which begin from a random config- 
uration and move around the discrete circle Z/iVZ following these dynamics: each 
particle waits an exponential amount of time to jump, it then chooses one of its near- 
est neighbours with equal probability. The exponential rate for the first particle to 
leave a particular site is a function c(-) of the number of particles at this site, and it 
is for this reason that the system is called the zero range model. Configurations of 
the system are denoted by 77; if we are currently at site x G Z/iVZ then the number 
of particles at the site is denoted by r](x). Next, we differentiate between particles 
by assigning to each one of k colours. The dynamics for each particle are the same 
as previously, however, in the colour process we keep track of the number of particles 
of each colour at all sites x of the discrete circle, and we denote this as r]i(x). If 
we ignore the colour of each particle then we obtain the colour-blind model, which 
simply keeps track of the total number of particles at each site x. We will say that 
the colour version contracts to the colour-blind process. 

The dynamics we have described conserve the total number of particles of each 
colour. Furthermore, the configurations have a family of invariant measures, indexed 
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by the average density vector p. These are usually referred to as the grand canonical 
measures, and we denote them here by p p . For the colour-blind model these contract 
to the measures \i p indexed by the total particle density. 

The above construction creates the model on a microscopic level. We, however, are 
interested in the system from a macroscopic viewpoint. In order to achieve this we 
re-scale space by 4. We then need to speed up time by iV 2 so as to obtain a non- 
trivial system evolution. This is the standard diffusive space-time re-scaling. We will 
denote the configurations of the re-scaled system by r] N and rjf for the colour-blind 
and /c-colour systems respectively. After re-scaling, the particles move around the 
discrete subset = {0, . . . , (N — l)/N} of the unit circle T. 

We first wish to understand the evolution of the macroscopic densities for the colour 
process. To do this, we study the asymptotic behaviour of the empirical densities 

We assume that the system is started in such a way so that the initial empirical 
densities correspond to some fixed macroscopic densities 

liniAr^oo n^ = Po(x)dx, (1.1) 

in the sense that this limit exists weakly on T. 

In the colour-blind model it is well known that (see for example |KLj ) that the 
macroscopic densities pt{x) evolve according to the equation 

dtp = \VD{p) Vp, p t \t=o = Po (1.2) 

with Po = Po given above. We call D(p) the bulk diffusion coefficient. Results of 
this type are known as hydrodynamic scaling limits in the literature. 

It is also known that the empirical measures Ylf t satisfy a law of large numbers 
and converge weakly to a non-random limit p\(x)dx (see CULL for a discussion). 
The vector p t , describing the limit, is the unique solution of the coupled non-linear 
differential equation 

dtp = \vD h {p) Vp, (1.3) 

with boundary condition p = {po}i=i- Here -Dfc(p) is the colour diffusion coefficient 
and is given by the formula 

(D-S)(p) 



D k (p) = S(p)I + 



P 



P 1 ••• P 1 



p k ... p k 



;i.4) 



where I is the k x k identity matrix and S(p) is the self diffusion coefficient which 
we now describe. 

Consider the zero range particle system where we have singled out, or tagged, one of 
the particles. The remaining particles are combined together to define a random envi- 
ronment of the tagged particle. This new system is a Markov process with identifiable 
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invariant measures. Using the methods of [KVj it is straightforward to show that this 
tagged particle started in equilibrium will converge to a diffusion with generator 

C = ±S(p)A. 

S(p) is called the self diffusion coefficient. One of the nice features of the zero 
range model is that we can calculate S(p) explicitly. Additionally, in the zero range 
case, we may express the difference between the bulk and self diffusion coefficients, 
D — S, as S'(p)p. We calculate these coefficients explicitly for the zero range model 
in Remark 13.11 We also use the notation x(p) f° denote the static compressibility. 
For the zero range model this is the same as the variance of the variable rj(x) under 
the invariant measure p p , and the following identity links the static compressibility 
with the bulk and self diffusion coefficients: 

S(p)p = X (p)-D(p). 

It is important to note that this is a particular property which holds for the zero 
range process, but not for general systems. 

The hydrodynamic scaling limit for the colour densities of an interacting particle 
system was first studied in | Qua for the symmetric simple exclusion process. This is a 
more difficult problem as the simple exclusion colour process is of non-gradient type. 

In this paper we study the central limit theorem associated to the law of large 
numbers described above. That is, we consider the fields defined by 

< Y»,f >= Vn\^ £ /(*)[<!(*)-**)]}■ (1-5) 

The main result is the following theorem, which identifies the asymptotic behaviour 
of Y t N . Let p t be the solution of (jl.3j) . Denote by P/v the measure on D{[0,T], H* 4 } 
induced by the stochastic field Y N . 

Theorem 1.1. Under the assumptions given in Section 3.3, P^ converges weakly in 
D{[0,T], H^_ 4 } to a coupled generalized Ornstein-Uhlenbeck process characterised by 
the coupled generalized stochastic differential equation 

dY t = ^{AD k (p t )}Y t dt + {VAl /2 (p t )}dW t (1.6) 

with initial condition Yq ~ Pq. In the above D k is the colour diffusion coefficient 
described previously, A k is given by 

A,(p)=™& ? {/ 1 ... )ft } ) (1.7) 

where x(p) ^ s the static compressibility, and D is again the bulk diffusion coefficient. 
W t is the k-dimensional Gaussian random field defined as a vector composed of k 
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independent copies of the Gaussian random field W[, with covariance 

E[(Wif){Wi,g)] = mm{s,t)<f,g>, 

for i — 1, . . . , k. 

The emergence of the colour diffusion coefficient in (jl.fij) is not surprising, as the 
equation may be obtained by formally linearizing the diffusion equation satisfied by 
the macroscopic density profiles. This notion is made precise in a result called the 
Boltzmann-Gibbs principle. 

Equation (jl.6|) contracts to the colour-blind density fluctuation field which we may 
write in the form 

dY t = l -/\D{ Pt )Y t dt + Vs/\xD]{p t )dW t . (1.8) 

We next consider the deviations of the colour fields away from Y t . To this end, define 
XJ\ = Y t l — pY t . The generalized SDE satisfied by U t may now be written as 

dU t = ^AS( Pt )U t dt + V^/S{pt)dZ t} (1.9) 

where Z t is a Gaussian process independent of Wt, such that 1 • Z t = 0. In terms of 
the original formulation we may express dZ\ as \fp~[dWl — Y^i=i V^^Wf, which 

is independent of Yli=i \f$,dWl = \fPtdW t by direct calculation of the covariances. 
Thus we have that Z t has covariance structure given by 

nrrmx{s,t) r f A / \ >. 

E[{Zl,f){Z>,g)] = I J f{x)g{x)pl{x)^-^ \dxdu. (1.10) 

Given (ll.9|) and (jl.Hjl . one may recover p.fijl . 

It is our belief that for general interacting particle systems, such as symmetric 
simple exclusion, the colour density fluctuations will be described by equations (|1.8|) 
and (|1.9|) . It is only the special relation [xD]p = S(p)p which allows the formula (|1.6|) 
to emerge for the zero range process 

Using equations ()1.6|) and (|1.8p . we may also write the colour density fields as 

* * Pt y Pt 

for i = 1 . . . k. The advantage of this formula is that we can see explicitly how the 
density fluctuation field for any colour interacts with the environment created by the 
colour-blind process. If our system was formed instead by non-interacting random 
walks the middle term in the above formula would be equal to zero. 

In the literature, results on the fluctuations of the hydrodynamic scaling limit are 
known as n on- equilibrium density fluctuations. When the system is started in an 
equilibrium measure, the terminology equilibrium density fluctuations is used. In 
fact, the term "density" is often omitted in the discussion; however, we introduce it 
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to emphasize the difference between these central limit theorems and fluctuations for 
the tagged particles. 

Next consider the empirical measure for the particle paths themselves. We de- 
note the paths of the diffusively re-scaled particles as X^(-), each taking values in 
-D([0, T], T). For each N we have n particles and we assume that 

11 

lim — = p (1.11) 

exists. The n particles are initially randomly distributed on the unit circle. Define 
the particle path empirical measure as 

Hn= J 5Xl (1-12) 

Il n is a random variable taking values in Mi{D([Q, T], T)}, the space of probability 
measures on D([0, T], T) endowed with the weak topology. That is, a realization of 
IT n is a measure obtained by assigning mass - to each observed particle path. 
A known result is that U n satisfies a law of large numbers 

n n -> n, (i.i3) 

where the limit II is a non-random element in the space Mi{D([0, T], T)}. In fact, II 
is the measure of a diffusion with generator 

\s( Pt )A, (1.14) 

where p t is the previously discussed macroscopic density and S is the self diffusion 
coefficient. The law of large numbers in (|1.13|) implies that a randomly chosen tagged 
particle has the asymptotic distribution II. It also implies that the particle becomes 
independent of the environment created by the system. This type of law of large 
numbers is often called propagation of chaos (see, for example, [Szn3j). Propagation 
of chaos implies a central limit theorem for the position of a randomly chosen tagged 
particle. 

Propagation of chaos was proved for the symmetric simple exclusion process in 
IQuat lR"ez] . First, the unique limit of Il n is identified using the colour hydrodynamic 
scaling limit Qua| . The fluctuation results of this paper imply that the hydrodynamic 



scaling limit for the zero range process satisfies equation (jl.3j) . Tightness in the space 
iWi{.D([0, T], T)} would complete the result. This last fact was proved in |Rezj for 
symmetric exclusion, and the method presented there may be applied to the zero 
range process. Large deviations associated to the above law of large numbers were 
also studied for the symmetric simple exclusion process in d > 3 QRV . 



The next goal is to study the associated central limit theorems. That is, we wish 
to understand the limiting behaviour of 

r n = v^(n n -n). 

In the case of independent random walks it is easy to see that the limiting distribution 
of T n is a mean zero Gaussian random field with the identity variance operator. When 
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there is dependence between the particles the limit remains Gaussian, but we expect 
an additional source of variation caused by the interaction. In Section |2] we discuss 
how one may obtain the above limit for the zero range process as a direct consequence 
of our non-equilibrium colour density fluctuation result. The type of convergence 
is convergence in finite dimensional distributions. We also calculate explicitly the 
variance for a class of test functions, and provide a formula for the inverse of the 
variance in the general case. 

Hydrodynamic scaling limits and associated fluctuations are both subjects of much 
interest in the literature. Scaling limits and equilibrium fluctuations are well under- 
stood for many models (see for example |KL| IML] IRos| \Uu\ IGKLj ) . However, there is 
still much work to be done to fully understand non-equilibrium fluctuations. To our 
knowledge there are currently only partial results for gradient models, and no known 
results in the non-gradient case. We adopt here Chang and Yau's |CYj proof for the 
Ginzburg-Landau model in dimension one to our zero range models. The proof relies 
on knowing the logarithmic Sobolev inequality for the inhomogeneous zero range pro- 
cess jjJanj . We apply this inequality in the proof of /c-colour non-equilibrium density 
fluctuations. The main idea behind the argument is that we may think of the zero 
range process for the i-ih colour as a process in a random environment, and our use 
of the logarithmic Sobolev inequality for the inhomogeneous setting reflects this no- 
tion. In order to make use of this result we make an additional technical assumption. 
For example, we could assume that the particle jump rate, c(k), is linear for large 
k (see Section 13.31 for a complete discussion of the assumptions) . We do not need 
this assumption elsewhere in this work, however, the full result holds only in this 
setting. The remainder of the work holds under the usual assumptions on the zero 
range process, namely 

(LG) swp k \c{k+l)-c(k)\ < oo (1.15) 
as well as a weak monotonicity condition 

(M) inf fc {c(A; + k ) - c(k)} > 0, (1.16) 

for some integer ko- 

Propagation of chaos and the associated fluctuations have been studied previously 
for systems with mean-field interactions, Sznl[ Szn2 IST| IAB| ITanj . However, the 



interactions between particles in the mean-field setting are weaker than for the zero- 
range process. The method presented here to obtain the result for finite dimensional 
distributions has not been previously applied to show a fluctuation result. It is based 
on the ideas developed in |Quat lR"ezl |QRV| . As the colour version of the symmetric 
simple exclusion process is of non-gradient type, at this time we cannot prove tagged 
particle fluctuations for this model. 

The outline of this paper is as follows. We begin with a discussion of the relationship 
between nonequilibrium density fluctuations for the colour version of the process and 
the central limit theorem for the particle paths. This is Section where we also 
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provide the formulae discussed above. In Sections El and IH we define notation, give 
a complete summary of the assumptions made, and give some preliminary results 
which we will use in the proof of the main result. Since the proof of Theorem 11.11 is 
quite involved, we first explain the result for the colour-blind model. This is done in 
Section In Section |H] we explain how to extend this to the colour version of the 
model. 

2. Fluctuations of Particle Path Empirical Measures 

A classical result from probability theory identifies the fluctuations of the empirical 
measure of independent random variables as a Gaussian field with identity covariance 
operator. More precisely, let Xi,X2,... be a collection of independent and identi- 
cally distributed random variables with values in some Polish space X and common 
distribution tt. We define the empirical measure, 7V n , as the random measure created 
by putting mass - at each of the n observed X^ random variables: 

The strong law of large numbers implies that 

in probability in the weak topology of probability measures on X. To study the 
fluctuations of this convergence we re-scale the quantity n n — n by ^pa. That is, we 
now consider the quantity 

Q n = ^/n(jj n - n). 
A consequence of the classical central limit theorem is that 

G n =^ G, 

in the sense of finite dimensional distributions, where G is the Gaussian random 
field with mean zero and identity covariance operator. That is, for each measurable 
function / : X — > ffi. in L 2 (tt) such that J fdir = we have 

< F n ,f > => <F,f>, 

where < T, / > is a real- valued Gaussian random variable with mean zero and vari- 
ance / pdu. 

The colour density fluctuations imply a similar result for the fluctuations of the 
empirical measure where the random variables are taken to be paths of the nearest 
neighbour zero-range particles after diffusive re-scaling. In particular, we have a sys- 
tem of n interacting particles with trajectories taking random values in D([0,T],T). 
One would expect that these fluctuations also converge to a Gaussian random field; 
however, the random variables under study are no longer independent and this non- 
trivial dependence structure would introduce an additional correlation in the limiting 
covariance. 



8 



HANNA K. JANKOWSKI 



For each N we have n particles and we assume p.lljl . The n particles are initially 
distributed on T and we assume that this initial distribution satisfies both a law of 
large numbers and a central limit theorem. That is, we assume that 

lim^oo i Yh=x ^xf (o) = po(x)dx 

exists weakly on T. Notice that this implies J po(x)dx = 1. We also assume that 

lim^oo (i Y2=\ S Xi(o) ~ po(x)dx) = T 

exists, where T is a random element of 7i_4. These assumptions are implied by 
assumptions (Dl) and (F2) of Section EHfl 

Define the empirical measure as in ()1.12j) . We next describe a heuristic approach 
used to identify its limit II. Consider the time-marginals of the process IT n . Define 

From Theorem 14.61 we have that H n (t) satisfies a law of large numbers and converges 
weakly to the solution of (|1.2|) . We would expect that a tagged particle started out of 
equilibrium would also converge to a diffusion with the self diffusion coefficient S(p t ) 
and a drift term caused by the evolution of the system towards equilibrium. That is, 
the generator of this diffusion should be 

C=~S{fH)A + b-V. (2.1) 

We next identify the drift b by noting that under the tagged limit we would still 
expect the time marginals of II to evolve according to the hydrodynamic diffusion 
equation (jl.2j) . That is, 

C*pt = • D( Pt )V Pt . 
Equating with the forward equation from the generator, we obtain 

b = ^-{VS(p t ) Pt - VD{p t )Vp t }. 

In the zero range process we have b = 0, as can be shown through direct computation 
of the above quantities. That is, for the zero range process, the limiting distribution 
II is equal to the measure P, where P is the law concentrated on C([0,T],T) of a 
drift-free diffusion process with generator given by 

C = \s( Pt )A. 



We wish to study the fluctuations of the previously described law of large numbers. 
That is, we center and re-scale as per usual, and study the quantity 

r n = v ^(n n -n). (2.2) 
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We have already identified IT as the measure P defined above. Theorem 11.11 implies 
that < r n , F > converges to a Gaussian random variable for smooth functions F : 
C([0, T], T) !-»• R of the form 

F(z(-)) = f(z(t 1 ),...,z(t k )), 

where t\ < . . . < tfc. Functions of this form correspond to the question "what are the 
particles doing at {ti, . . . We answer the question by colouring the particles 

according to their behaviour at these times, and hence, we may use the fluctuation 
results for the colour densities. 

To describe this consider first the simplest case of 

F«-)) = /(*)• 

Denote by Vf s the semi-group associated to the nonpositive operator ^D(p t )A. The 
results of Section El tell us that for this choice of F, < T n , F > converges to the 
random variable described by 

E[e^ F >] = MV°J)eM~E[{j\Kf,dZ s )\ ]}, (2.3) 

where 7?o(/) is the characteristic function of the random variable < T , / > . In the 
above Z t denotes the generalized Gaussian process with covariance 

rmin(s,t) 

E P [(Z t ,f)(Z s ,g)]= / <Vf,[xD]{Pu)Vg>du. (2.4) 

Jo 

Indeed, this is simply a restatement of Theorem 15.81 
Consider next the case of 

F(z(-)) = Ilzih) E B 1 ]f(t), (2.5) 

where t>t\ and E>\ is an open subset of T. Following [Rez] , we define 



Pl(x) = p tl {x) - pl^x). 



(2.6) 



Let Vf s denote the semi-group associated to the equation d t h = \S(p t )/\h. We also 



define A\ as the operator acting on functions h as - — — pi^A/i. For the zero-range 
process this becomes S"(p)p 1 A. Lastly, as in the above, we define Z\ to be two 
independent generalized Gaussian processes with covariance 

/■min(s,t) i 

E P [(Zi,h)(Zt,g)} = < V/, [ x D}( Pu )^Vg > du. 

JO Pu 
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The results of Section |H1 imply that for the choice of F given as in (|2.5|) . < r n , F > 
converges weakly to the random variable described by 



(2.7) 



r 1 ^ 

X exp{--E 



(P t s J,dzl) + f\vl u i Bl vl t j,dzl + dzl) 

Jo 

+ f f^AlVfj^dZl + dZDds 

Jti Jo 



}• 



Notice that in the above formula the distribution of Z\ + Z^ is the same as the 
distribution of the field Z t with covariance given in (|2.4j) . We chose to leave the for- 
mula in this format to emphasize the relationship between the terms. This formula 
is obtained by formally solving the generalized Ornstein-Uhlenbeck stochastic differ- 
ential equation. That we may do this is made rigorous by Theorem 15 . 81 together with 
Theorem 16.71 

We may now repeat the above to obtain convergence for general functions F for 
any number of k. That is, 

F(z(-)) = UlMU) G B t ]f{z(t)). (2.8) 

The ti are fixed and increasing times in [0, T] with ty. < t < T, By. are fixed measurable 
subsets of T, and / is a smooth function. Notice that each step doubles the number 
of colours used, and we have 2 k colours for the general case. 

Except for the initial measure 7r of r , our limits are Gaussian. If we assume that 
7To is also a Gaussian measure on 7i_4 then we would have that V is truly a Gaussian 
field. We do not make this assumption; however, from this point on we will say that 
T is a Gaussian random field and omit the discrepancy caused by the initial measure. 

It follows from the above that we have proved weak convergence of 

<r„,F>^<r,F> 

in the sense of finite dimensional distributions for a large class of functions on 
C([0, T],T). We know that T is a mean zero Gaussian random field. We would 
like to identify explicitly the covariance operator of V for a general function G. 
For a field G : C([0, T], T) — > R define 

q(t,x) = E P [G(X(-))\X t = x}p t (x), 

where p% is again the solution to ([1.2)1 . Let 7* denote the quantity 

7t ' j " lx-D]( Pt (x)) ' 

where Q is the functional 

g(t,p,q)=d t q-^A[D(p)q]. 
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Lastly, define the operator A*G = J Q Y(t, X t )dX t . 

Theorem 2.1. For any mean-zero function G in C& [C([0,T], T),R] into the real 
numbers, < T,G > is a real-valued Gaussian random variable described by 

E[e t<r ' G> ] = vr (Fo) e ~^(G) e -iQ s (G)_ (2 g) 

7To is the characteristic function of the initial field T and F = Ep[G(X(-))\X ] is 
the projection of the function G onto the initial field. Qd(G) = Ep[G Q~oG] may be 
identified through the operators 0£> on a test function F corresponding to the operator 
on a field G given by the quadratic form 

E P [(A*G) 2 ]. (2.10) 

Similarly, for Qs(G) = E P [G Q^G] may be identified as the inverse of the quadratic 
form 

E P [{G - A*G) 2 ]. (2.11) 

Remark 2.2. In the case when we are working with a system of independent random 
walks (that is, we choose c(k) = k), the above formula becomes Ep[G 2 } as expected. 

Proof. It is enough to establish equivalence with the limiting behaviour of < T, F > 
for functions F of the form in ()2.8|) . 

For any function G, the behaviour of < T, G > may be split into three components. 
First, we have the behaviour of < T, Ep[G(-)\Xq] >, which is independent and given 
by the initial field measure ir . Next, we consider the evolution of the marginals 

<T,Ep[G{-)\X t }- Ep[G(-)\X ] >, 

whose variance is described by the operator A*G = A*[Ep[G(-)\X t ] — Ep[G(-)\X }}. 
This is the same as (jl.8|) and by formula (|2.9|) is independent of what remains. 
The remaining portion is the behaviour of < T,G > less the time evolution of 
< T,E P [G(-)\X t ] -E P [G(-)\X ] > . By (Q and IjgTQ this is Gaussian with identity 
variance operator under the measure P. For functions of the form given in (J2.8)) . 
this is simply (jl.9j) . which is independent of (|1.8|) . That is, consider the case of 
F(z(-)) = f(zT)^-B(z s ). Using ()2.6|) and the propagation of chaos results of jRezj . we 
have that E P [f(X T )I B (X s )\X t = x]p t {x) = E P [f(X T )\X t = x]p l t {x), as in the latter 
case the indicator function is simply equal to one. It thus remains to argue that the 

appropriate density in this case is given by p\{x) jl — j to match (|1.10|) . This 

follows from noting that Pt(%) = Ym=i p\{x)l\paxtic\e is of type i], for any number of 
colours. Thus, we may identify G^ 1 via the variance given in (|2.3jl . and Q^ 1 is simply 
the identity on G - A*G. □ 

Theorem 12.11 tells us that the limiting fluctuations of propagation of chaos for our 
model separate into the bulk density fluctuations and the independent fluctuations of 
the remaining deviations. Fluctuations of this type have been studied previously for 
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systems with mean-field interactions in |Sznl[ IHzn2[ ISTj IAB[ ITanj . where the limiting 

behaviour is quite different. For these models, the variance operator may be written 
in the form (I — A) 2 , with A described explicitly in terms of Malliavin derivatives of 
P. The physical source of this form for the variance operator is most easily seen in the 
coupling approach developed in |Sznlj . We also note that the convergence obtained 
in |Sznl| ISzn2| IST| IAE>1 ITanj for the fluctuation results is weak convergence in finite 
dimensional distributions only. 



3. Notation and Assumptions. 

3.1. Notation. We denote the set {0, 1/N, . . . , (JV - 1)/N} as T N , with addition 
defined as addition modulo one, and the continuous unit circle as T. 

Given a subset A of T N or of Z/7VZ we write AV x£ \f (x) to denote the average of 
the function / inside the box A, that is, 

AV xeA f(x) = 

For a metric space X, C P (X) stands for the space of real-valued functions on X with 
p continuous derivatives, and Cb(X) the space of bounded and continuous functions 
on X. 

For z in Z define h z : T i— > R by e z (x) = \/2cos(2irzu) for z positive, e z (x) = 
y/2sia(27rzu) for z negative, and eo(x) = 1. The collection {e z } is an orthonormal 
basis of L 2 (T). < •, • > denotes the inner product of L 2 (T). The functions e z are also 
eigenfunctions of the operator I — A , with eigenvalues 7 Z = 1 + 4ttz 2 . We use the 
notation Ti m to denote the Hilbert space formed by taking the completion of C°°(T), 
the space of infinitely different iable real functions on T, under the inner product 

<f,9> m =< f, (I-A) m g>. 

For each positive integer k we denote by H- m the dual of Ti m relative to the inner 
product < •, • >. 

We define the k-fold product of these space to be 

T~£_ m = T~t-m x 7~L-m x ... x T~t_ m . 

That is, for / G 7i^ m we have / = {fi, ■ ■ ■ , fk} and each fi is an element of TC- m - 
The norm ||/||_ TO| fe is defined as 

llAll-m + II/2II-TO + • • • + H/fcll-m- 

We also define the norm 1 1 / 1 1 * for a function defined on the discrete unit circle and 
satisfying the condition J2 x &t n f( x ) = ® as 

11/11. = ^ J2 f( X K- A Nf\(^ 
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where A at denotes the discrete Laplacian 

[A N h](x) = N 2 {h(x + 1) - 2h(x) + h{x - 1)}. 

We denote by < •, • >* the associated inner product. The traditional notation for this 
norm is also || • ||_i, which we do not use to avoid confusion with the Sobolev spaces 
defined above. We also define the product norm to be 

1 1/| \*,k — H/ill* + 1 1/2 1 1* + • • • + 1 1 /fell*) 

for / = {/i,...,A}. 

3.2. The Model. The class of zero range particles we describe is the symmetric 
nearest neighbour zero range interacting particle system in dimension one. 

The Colour-Blind Process. First consider the evolution of the number of particles 
at each site. If a particle moves from site x to site y the configuration rj changes to 
7] x ' y where 

{rj(z) — 1 if z = x, 
rj(z) + l if z = y, 
7]{z) otherwise. 

It does this at rate ^c(rj(x)). The system is a Markov process and we can write down 
its generator L as 

( L f)(v) = \Y< c ^{x))[f{ri x > y ) - M\ (3.1) 

where x ~ y denotes nearest neighbours of Z/iVZ. 

For ip > define the partition function Z((p) = J2k>o "my." Under the assumptions 
(LG) and (M), the radius of convergence for Z(-) is infinite. Fix < (p < 00 and 
denote by fi^ the product measure on X with marginals 

» Mx) = k) = z$4W 

where c(k)\ = rii<m<fc c ( m ) if ^ > and c(0)! = 1. The family of measures fi v 
indexed by <p are stationary and reversible for the zero-range process. If we define the 
product measures on the infinite lattice Z, and not on T, these measures represent 
the full set of extremal invariant measures for the system |Andj . 



The measures fi^, are referred to as the grand canonical measures in the physics 
literature. They are not ergodic, however, as we have already pointed out that the 
system evolution preserves the total number of particles. If we condition on the aver- 
age number of particles, y, we do obtain stationary, reversible, and ergodic measures: 
Hpivlv = y)- Equivalently, we may also condition on the total number of particles. 
These conditional measures are called the canonical ensembles. 

To emphasize that these measures are defined on the discrete unit circle of size N 
we will use the notation v N ^ y . That is, v Ny = ^(77)77 = y). If we are considering the 
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measures for configurations restricted to a subset of the circle, and this subset is of 
size K, we will use the notation v K y . Because of the homogeneity of the system there 
is no ambiguity in the notation. 

Let p(ip) = En v [r)(x)] denote the average density of particles. By straightforward 
computation we obtain the following identities 

p(<p) = V Z\ V )Z-\ V ), (3.2) 
f/(<p) = V - 1 E flv [(r ] (x)-p( V )) 2 }>0. (3.3) 

As the variance must be strictly positive for ip > 0, we conclude that p(ip) is a strictly 
increasing function and hence invertible. Because p has the natural interpretation of 
the density, we will fix p and think of ip as <f(p). We shall then index the invariant 
measure by p: p p . Notice also that E Pp [c(rj(x))] = y?(p). 

The static compressibility mentioned in the introduction is in general defined as 

Y. x ei/Ni E ^Hx)]r](Q)l 

where we use the notation E[f; g] to denote the covariance of the functions / and g. 

We are now in the position where we may express the self and bulk diffusion 
coefficients, as well as the static compressibility explicitly. 

Remark 3.1. For the zero range process we have: 

S(p) = D(p) = <p'(p), and X (p) = ^f\. 

p tp'{p) 

This follows from the above discussion of the invariant measure p p . 

We next define the Dirichlet form, D flp (f) = E Pp [f(—L)f]. Since p p is reversible 
for the dynamics, this is equivalent to 

x~y 

The same identity holds when the expectation is taken with respect to the canonical 
ensembles, in which case we denote the Dirichlet form as D UN {f). We may also 
restrict the Dirichlet form to a subset of the discrete circle, Ak, where |A| = K. Here 
we again use the notation D UKy (f). 

D » K M) = l E E ^Mv(x))[f(v x ' y ) - mn (3.5) 

The Colour Process. To simplify notation we define the model for the case when 
there are only two colours. The extension to general k is immediate. 

Imagine that the zero range process is made up of two different colours of particles. 
The two types of particles are mechanically identical to the regular zero range process 
particles, but the two-colour process keeps track of the two types of particles as the 
system evolves. That is, we are now studying the evolution of the number of particles 
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of each type at site x in Z/iVZ. Its elements will typically be denoted by the pair of 
configurations 77= (771,772) £ N z / ra x N z/iVZ . The dynamics for each particle are the 
same as in the zero range process. Thus, the first particle of colour i to jump from 
site x does so at rate 

_ 7ft (z) 0(771(3;))+ 772(3;)) _ , , c(77(x)) 
C ^ X >>- ( Vl (x)+ V2 (x)) - m{X) V (x) ' 

where r)(x) = rji(x) + 772 (x). We write down the infinitesimal generator for the two 
colour process 

(£/)(*?) = \i2Y,«W x ))\fW v ))-fto)i ( 3 - 6 ) 

i=i 

Here r)*' y denotes the configuration obtained from 77 by moving one particle of colour 
% from site x to site y. Note that if the function / is "blind" to the particle colour, i.e. 
fiv) = f(Vi + then Lf is equivalent to the generator for the previously defined 
zero-range process. Because of this contraction, we shall not use a different notation 
for the generators of the two processes. 

We now define the grand canonical measures and the canonical ensembles for this 
process. Fix cp = {cpi,cp 2 } > and denote by \x v the product measure on X 2 with 
marginals 

1 ( \ u ( \ \ V 9 ! ^2 fk + m\ 1 

/MW) = k : rj 2 {x) =m) = — — — -. 

c(k + m)! \ k J Z(cp 1 + ip 2 ) 

Z(cp) is the partition function defined previously. We also have that 

AVM^) = k\r] 2 (x) = m) = x (3.7) 

Cm{k)\Z m [ifi) 

where c m (k) = kc ^+m^ an d Z m (cpi) is the associated partition function. The family 
of measures indexed by cp are stationary and reversible for the two-colour zero range 
process. 

Let p l {<p>) = [?7j(a;)] denote the density of particles of the i th colour. Notice that 

p = ~zW and 

d, Pj p i = <Pj 1 E llv [ri i (x)ri j (x)-p i p']. 

In particular, this implies that the Jacobian of the transformation cp 1— > p = {p^p 2 } 
has determinant strictly positive for cp > 0. As before, we choose to index the 
invariant measure /i^ by the pair (p 1 , p 2 ), (that is, work with p = p p i p 2), where we now 
consider cp { = cpi(p). The function ipi can be recovered through cp { = E fMv> [c i ( y r](x)] = 

P % ^~v) = Since c(r}i(x) + r] 2 (x)) = 01(77(2;)) + c 2 (r)(x)) : we obtain that cp = 

<fi + cp 2 — E^ v [0(771(2;) + 772(2;))] for the colour-blind model. 



16 



HANNA K. JANKOWSKI 



As before we define the canonical ensembles to be the measures // v conditioned on 
the average density of both colours of particles, which we denote as y — {2/1,1/2}, 

VN, y (r}) = vMVi = Vi,i= 1, 2). 

If we consider this measure only for configurations restricted to a subset of the circle, 
where the subset is of size K, we will use the notation VK,y 

The Dirichlet form for the two-colour version, (/) = E^ [/(— L)f], is then 

lY.^M'nmm*) - /N] 2 ] +\Y^EM^))[m v ) - m\*\ 

x~y x~y 

= Dl(f)+Dl p (f). (3.8) 

The same identity holds for the canonical ensembles, in which case we have the 
Dirichlet form 

D, N Jf) = E m Jf(-L)f] = Dl N Jf) + Dl N Jf). 

We will use the notation D„ K (f) as in the previous section, when we have restricted 
the Dirichlet form to a subset of size K. 

We denote the expectation with respect to the measure induced by the zero range 
process started in equilibrium as E E q and as E NE q if the process is started out of 
equilibrium. Because of the previously noted contraction, there is no contradiction 
in using this notation for both the colour and colour-blind models. 

3.3. Assumptions. We separate the assumption into several subsections. 
On the Rate Function. 

We make the assume that the rate function c(-) satisfies (LG) and (M) of (|1.15)) and 
Assumption (LG) is necessary to ensure that the zero range process is well 
defined on the infinite lattice |Andj . Condition (M) rules out the cases, such as the 
queueing system corresponding to c(k) = I(k > 1), where the spectral gap is known 
to depend on the density of particles. 

A key ingredient in the proof of nonequilibrium density fluctuations is the logarith- 
mic Sobolev inequality. In order to make use of this tool in the inhomogeneous zero 
range process we make an additional assumption. We make the assumption only to 
be able to quote this result, and it is not required in the remainder of the work. 

(E). Recall the conditional grand canonical measure from (|3.7|) . We assume that for 
every N e N there exists a finite positive constant C = C(Nq) such that 



m) <C, (3.9) 



jx) = r 

, x-eA 

holds uniformly over all r G N and configurations rj on a subset A C Z such that 
|A| = N . 
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Because of the contraction principle, it is sufficient to make this assumption for 
only two colours. Also, by symmetry, it follows that this assumption also holds for 
any colour, conditioning on the configuration of the remaining particles. This is a 
very technical condition, and it is satisfied if instead we make, for example, one of 
the following assumptions. 

(El). There exists a large constant K , and a positive constant 9 such that for all 

k > Kq the rate function c satisfies c(k) = 9k. 
(E2). There exists a large constant K , and two positive constants 9i and 9i such 

that for all k > K the rate function c satisfies for all x 

(t.\ — / ^ ^ ^ ^ is odd, 
1 #2 k if A: is even. 

Naturally, many other variations on these exist. For more details on this assumption 
see jJanj. 

Remark 3.2. The assumptions (LG) and (M) imply that there exist finite positive 
constants c\ and c<i such that for all k c\k < c{k) < c^k. 



On Initial Density. 

Let Co denote a sequence of probability measures on the k-fold product of N Tn . For 
each fixed N, £q stands for the initial measure of the colour zero range process rjf 
with k colours. 

(Dl). We assume that for each fixed k the are associated to some fixed density 
profile po : T fc — >• [0, 1] (where po = {po}«=i,. -,fc ) i n ^ e sense that 



lim E, 



Mi 



i fWN)n?{*) - [ 



f{x)p l {x)dx) 



for any continuous function / G C(T) and for each i from 1 to k . 

(D2). We also make the assumption that 

lim sup [AV x {rj N (x)} 2 ] < C. 

N^oo 

On Initial Entropy. 

(HI). We assume that the relative entropy of the intitial measure grows at most 
linearly with respect to an invariant measure. That is, for any k, assume that there 
exists cp* > and a constant C = C(k) such that 

H(CM<CN. 

Remark 3.3. Using the entropy inequality we can show that if the above bound is 
satisfied for one vector cp* then it is also satisfied for each ip > 0. 
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On Initial Fluctuations. 

Let Yj^ = {Yqj_, . . . , Y^ k }, where each denotes the initial density fluctuation field 

Y Q N l (x) = VN{ V ?(x)-pl ) (x)}. 
Also let Y^ denote the centered version of Y^. That is 

Yol = JN{ri?{x) ~ Pl{x) - AV x {r^(x) - „•(*))}. (3.10) 

(Fl). We assume the following 

(a) . 

limsup^[7V-^||Y ^||y =0 

(b) . 

]imsap E C N[N 1/2 \rji - AV x£T Npi\ 2 } = 0, for all i = 1, . . . , A;. 

JV->oo 

(F2). Let M(H k L m ) denote the space of measures on H. k _ m . Assume that the initial 
fluctuation field has a weak limit in M(H k _ m ) for m = 4 . That is, if we denote by 
the measure induced on M(7i fc m ) by the initial field Y t N , then there exists a unique 
measure P such that =^ P in the weak topology of M(H k _ m ). 

4. Preliminary Results 
We discuss in this section some results used in the remainder of this work. 

4.1. Connection with Inhomogeneous Zero Range . The symmetric inhomo- 
geneous zero range process has dynamics which are the same as of the homogeneous 
process, except that the jump rate for a particle at site x now also depends on the 
site x. That is, denote the rates at site x as c x (-), and £(x) the number of particles 
currently at x. Then the first particle jumps from x at rate c x (£(x)), and it jumps 
to one of its nearest neighbours with equal probability. We assume that these rates 
satisfy the bounds (LG) and (M) uniformly in the environment. 

(LG*) sup \c x (k + 1) — c x (k)\ < ai < oo 
(AT) inf \c x (k + k ) - c x (k)\ > a 2 > 
for some constants di, a 2 , k . The generator of this process is 

(Lf)(v) = \ E cMx))[f(e' y )-f(0] (4-1) 
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The process has invariant grand canonical measures which are product measures with 
marginals 

where Z x (ip) is the appropriate partition function. The process is reversible with 
respect to these measures, as well as the canonical ensembles, obtained by conditioning 
on the average density of particles. We denote the canonical ensembles on a subset 

A as z>a i3/ (0 = £a,^(£IC = y). 

Notice that the mean density varies depending on site in the inhomogeneous setting. 
However, we may still consider (p as an invertible function of the overall density 
]X\Y, x€ a E va,Av(s)} for every fixed A. Similarly, we define a 2 x {ip) = £ Aa Jr^x); 7?(x)], 
and 

' ' xeA 

The Dirichlet form, with expectation taken with respect to the canonical ensembles 
on A K (x) = {y : \x — y\ < K}, shall be denoted as 

= \E»^ IE.,..., A( , : , : 'iii.r))i/iC'-) " /(0) 




We make the additional assumption that the inhomogeneous zero range satisfies 
the following property. For every Nq > 2, there exists a constant C = C(No) such 
that 

(E*) i < y/fHAM-fo) I > ^ =r)<C (4.2) 

uniformly in r G N, and sets A for any Nq 

Relationship with Colour Zero Range. Consider the colour zero range process 
where all except the particles of the first colour have been "frozen" . By the contraction 
principle, we may consider only the case of two colours without loss of generality. 

To this end, fix a configuration i]2 and let rji evolve as though it was a single-colour 
zero range process which uses the non-homogenous rate function 

This is a particular example of a non-homogeneous zero range process. The invariant 
measure for this process is clearly the product measure with marginals given in ([3.7)1 . 
the two-colour process invariant measure conditioned on the second colour configu- 
ration. It is not difficult to show that conditions (LG), (M) and (E) imply that the 
measures fl x (k) satisfy conditions (LG*), (M*) and(E*). 

The Dirichlet form (with respect to the canonical ensemble) in this case is simply 

d :L(a „M) = E^ ix) ,\ m) [E^a^ c*(Vi(x))(f(rii' y ) ~ f(v)) 2 ] ■ (4.3) 
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We shall use the results of [Jan on the inhomogeneous zero range process to state 
some useful facts about the conditional colour zero range. In what follows we shall 
state the results for colour 1 conditioning on colour 2, or equivalently, on all of the 
other colours. By symmetry, the results are valid for any colour. 



4.2. Moment Bounds. The following two lemmas are proved in [LSVJ. 

Lemma 4.1. There exist constants C\ and C 2 , which depend only on the values ai, a% 
and ko, such that 

< d < < C 2 < 00. 

P 

Lemma 4.2. For all k > 1, there exists a finite constant C(k) such that 

m 2k (p) < C(k) a 2k (p) 

for all p > p > 0, where rrij(p) denotes the j th moment ofrj(x) under the distribution 
p p . 

Next we prove a simple but useful property of the functions <p(p) and S(p) = "^y-'- 

Proposition 4.3 (Lipschitz properties). The following functions are Lipschitz: 

(i) Vifi) 

(ii) ifi(p) as a function of p % 
(hi) S(p) 

The first two functions are also strictly increasing. 

Proof. The proofs follow from direct calculations of the derivatives, as well as previ- 
ously computed bounds: 

(i) ^) = %■ 

(U) d p ^ l (p)=is(p)+^'( P ). 

(iii) First of all notice that if c(k) = k then S(p) = 1 and there is nothing to do. 
Otherwise, S'(p)p = ip'(p) — S(p) which implies that S'(p) is bounded below 
as long as p is bounded below. To finish we need to examine the limit of 
S'(p) as p — > 0. This follows from noting that 

limS"(p) = -lim^(p) 2 ^ x lim J W> W> W> W){ 

□ 



p-^o ' P ^o J a 2 P ^o [ Z 2 (ip) 



4.3. Spectral Inequalities. When the underlying measure p is one of the canonical 
ensembles Vx,y we will write the entropy as Hx, y {f) for ease of notation. The following 
inequality was proved by Dai Pra and Posta }DPPlj . |DPP2j under assumption (LG) 
and (M). 
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Theorem 4.4 (Logarithmic Sobolev Inequality). There exists a constant Cls = 
CLs{ a ii a 2iko) such that 

H K ,yU)<CLsK 2 D UK ^f) 

holds for any K, y and positive function f . 

If we make the additional assumption (E), in light of the discussion of section l4~T1 
and the results of [lanj, we obtain a logarithmic Sobolev bound for the conditioned 
zero range process as well. 

Theorem 4.5 (Logarithmic Sobolev Inequality for Conditioned Zero Range). Assume 
that conditions (LG), (M), and (E) hold. Then there exists a constant Cls — 
CLs{0'i, a '2,ko) such that 

HVWkMv*)) < C LS K 2 Eft^Vf) 

holds for any Ak(x), yi, density f on BJ Aif W ; and configuration of particles of the 
second colour i]2 . 

4.4. Hydrodynamic Scaling Limits. 

Theorem 4.6. Under assumptions (LG), (M), (Dl) and (D2) we have that for every 
t < T, for every continuous function g : T d — > M. and for every 5 > 0, 



lim P, 



N^oo 



i ^ g(x/N) Vt (x) - f g(6) Pt (6)d6 



> 5 







where pt{0) is the unique weak solution of the non-linear heat equation M.ty) . Here 
we use P(-n to denote the measure of the zero range process with initial measure . 

Remark 4.7. It is possible to prove this under reduced assumptions on the rate func- 
tion. It is sufficient that c(k) > 8k for all k > and some positive constant 9. 

The proof of the above fact appears in Kipnis and Landim ( |KLj ) and is based on the 
entropy method first developed by Guo, Papanicolaou, and Varadhan |GPVj . 

Remark 4.8. Theorem M . 1\ implies the hydrodynamic scaling limit result for the colour 
version of the zero range process. However, it is possible to prove this result indepen- 
dently and under fewer assumptions. Most notably, assumption (E) is not necessary. 
We may use, for example, the nonhomogeneous spectral gap from jjanj. 



4.5. Uniform Local Limit Theorems. For m > 0, denote by H m (x) the Hermite 
polynomial of degree m: 



H m {x) 



rf.1 Jm 

j iX- . Li/ 

exp( -T } ^ 
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Let go(x) denote the density of a standard normal random variable, and define for 

km 



9j(x) = g (x)^H j+2a (x) Yl 



1 



m=l 



k m \ \ {m + 2)\o m + 2 



(4.4) 



where the sum is taken over all nonnegative integer solutions {ki} J l=1 and a such that 
k\ + 2&2 + . . . + jkj = j and k\ + hi + . . . + kj = a. 

Theorem 4.9. For all p* > and J G N, there exist finite constants E Q = E (p*, J) 
and A = A(p*, J) such that 

J -2 



£>!( 

.x<=A 



X 



Np + aVNz 



3=0 



< 



E 



{a 2 (p)N)( J - 1 )/ 2 



uniformly over z and over all parameters p* > p > A/N. 

The above is an Edgeworth expansion for a lattice distribution, valid uniformly for 
the family of measures p p , the first as A/N < p < p . Its proof appears in |LSVj We 
also make use of the same result proved in |Janj for the inhomogeneous process. We 
state the result for the conditioned zero range process. In the theorem we have that 
\A\ = N. 

Theorem 4.10. For all p* > and J e N, there exist finite constants E = E (p*, J) 
and A = A(p*, J) such that 

J-2 



V2 



j=0 



< 



E 







(J-l)/2 



Y]rn(x) = Np 1 + aVNz 
_xeA 

uniformly over z, configurations of r]2, and over all parameters A/N < p 1 < p*. 

5. Density Fluctuations in Non-equilibrium 

Let p t be the unique solution to the partial differential equation (|1.2|) . Define the 
density fluctuation field for the colour-blind model as: 



<Y t N J>=VN{±J2f^( 



Pt(x)}}. 



(5.1) 



N 



Let Pjv denote the measure of on D{[0, T], H-4} induced by the stochastic field Y 
with initial measure . 

Theorem 5.1. Under the assumptions of Section \S.,% except for (E), Pn converges 
weakly in D{[0, T], i/_4} to a generalized Ornstein-Uhlenbeck process characterised by 
the generalized stochastic differential equation 



dY t = ^{A<p'( Pt )}Y t + {Vv^(a)}«. 



(5.2) 



NONEQUILIBRIUM COLOUR FLUCTUATIONS 



23 



where Wt is the Gaussian random field (Wt, f) with covariance 
COV[(W t , f) (W s , g)} = mm(s, t)<f,g>. 

The process Wt is also known as a generalized Brownian motion which is "Brown- 
ian" in time and "white" in space. 

Before proceeding with the proof of the above we make the quick remark that non- 
equilibrium density fluctuations have been proved for a specific class of zero range 
processes in |FPVj . The class of processes in their work is very special; the rate 
function they consider is c(k) = I[k > 1]. In this very unique case there exists a spe- 
cial equivalence between this zero range process and the symmetric simple exclusion 
process. In this setting, the martingale equations are closed in the field Y t N and no 
Boltzmann-Gibbs argument in necessary. 

5.1. Outline of Proof. The proof is divided into several main steps. We first identify 
the drift and quadratic variation of the limiting field Y t . We then show that there 
is only one measure which solves this martingale problem, and that the sequence of 
measures P N is tight. We first turn our attention to the drift and quadratic variation 
martingales. Suppose that P/v fe is a convergent subsequence of P/v, and for ease of 
notation, we denote it again by P/y. Let P denote its weak limit. 

Identifying the drift and quadratic variation of Y t which matches the equation (|5.2j) 
is equivalent to the statement that for any test function / 

< M t , f > = < Y t , f > - < Y , f > - < Y s , ~<p'(p s (-))Af > ds, (5.3) 
is a martingale under P with quadratic variation given by 

E P [(< M t ,f>-< M s J>f\T s \ = J*Jip( Pu (x))[Vf(x)} 2 dxdu. (5.4) 
We next compare these to the martingales under the measures P, 



AT- 



By direct calculation we know that M/ defined below is a martingale under P 1 



N 



for any N. 

< M t N , f >=< Y t N , f> - <Yo N J>-jt< Y S N , A N f > ds (5.5) 
where Y denotes the fluctuation field \/N \c(r] N (x)) — <p(p(x))). 



Comment on notation: There is a slight abuse of notation in the above as we use 
the same inner product notation < f,g > for functions /, g defined on T, where 
it is equal to J fgdx, and for functions /, g defined on where it is equal to 
7f J2xgt n fi. x )di. x )- We do this only to simplify the notation. 

A similar calculation gives the quadratic variation of < , f >: under P N , N^(f) 
defined below is a martingale. 



<M t N j> 2 ~ I a\;__^ 



[V N f(x)] 2 ds. (5.6) 
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It hence remains to show that the limits of and Nj* are consistent with (|5.Hj) 
and (|5.4j) . This is much easier for Nf*(f). The difficulty which arises in is that 
expression (|5.5|) is not closed in the field K^. This will prove to be the main obstacle 
to overcome in the proof. We start with the fact that < M t , f > is a martingale, 
and hence we know that for all N 

E Pn [< M t N , f>U]= E Pn [< Mf, / > U], (5.7) 

for all bounded, .^-measurable U. Here, T t denotes the cx-algebra on D{[0,T], H^} 
generated by F s (f) for s < t and / G C°°(T), where F G D{[0, T], #_ 4 }. We need to 
show 

E P [< M t , f>U}= E P [< M„ f>U], (5.8) 

Comparing ()5.3|) with ()5.5|) we see that we need to replace the field Y t N with ip'(p t )Y t N . 
This type of result is known as the Boltzmann-Gibbs principle in the literature. 

The remainder of this section will be divided in the following manner. We begin by 
showing that every weak limit of the measures Pjv solves the martingale problem. This 
is subdivided into identifying the asymptotic drift and then its quadratic variation. 
We first consider the drift, and begin with the proof of the Boltzmann-Gibbs principle. 
We handle the quadratic variation in Section 1531 In Section I5~4l we discuss the theory 
of Holley and Stroock with states that there is only one solution to the martingale 
problem. Section 1531 is dedicated to tightness of the measures P . These results in 
combination as outlined above prove Theorem 15.11 



5.2. Identifying the Drift Martingale. In light of the preceding discussion in the 
introduction, to show that under P the drift martingale is identified through (J5.3j) . it 
remains to prove the following. 



Theorem 5.2 (Boltzmann-Gibbs Principle.). For functions f G C^T] 



lim Eneq 



< /, VN(c( V ») - <p(p s ) - tp'(p s M ~ Ps)) > ds 



N 



0. (5.9) 



This section is dedicated to the proof of this result. 

The first step of the proof is to replace pt(x), the solution of (jl.2|) with the solution 
to the discretized version 



d t p = -A N (p(p), 



(5.10) 



with initial conditions Pt{x)\t=o = Po{x) for x in T N . The difference between the two 
solutions is of order and hence does not affect (|5.9j) . RMJ. Because of this fact 
and in order to simplify notation, we continue to denote the solution of f)5.10j) as p in 
the remainder of this section. 
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We begin by re-writing the field \/N{c{r]^) — <p(pt) — ^'(p^ivf ~ Pt)) as the sum 
of five separate parts: 

$f it) = Vn[t£(x) - c K {v?{x)) ~ ^{pt{x)){rg{x) - mf{x)}\ 

$»(t) = VN[c K (^(x))-^mf(x))]I[mf(x)<R] 

*»{t) = y/N[<p{m*(x)) - <p(pt{x)) - y' {p t {x)){m? {x) - p t (x)}]l[m? (x) < R] 

$f (t) = y/N[-ip(fH(x)) - V \p t {x)){mf{x) - p t (x)}]I[mf (x) > R] 

$f (t) = y/Nc K {x)l[mf{x) > R]. (5.11) 

Here, mf(x) = AV\ v - x \< K / N r^(y) and similarly, c K (x) = A V\ V „ X \< K/N c(r]^(y)). 
To prove Theorem ()5.2j) . we thus need to show that 

Eneq[\H <f,*?(s)> ds\]^0 (5.12) 

for for each i. We begin by proving several lemmas, which summarize the two main 
components of the Boltzmann-Gibbs principle. The argument is essentially a Taylor 
argument, and the local equilibrium principle gives us the first term in the expansion, 
while "equivalence of solutions" allows us to control the resulting error term. 

5.2.1. Local Equilibrium Principle. The first lemma we need proves that over boxes 
of microscopic size K we are able to replace the average of a function with the 
expectation of said function at the average density m K (x). We take K « N. It 
turns out in the end that the correct scaling for this result is to take K to be slightly 
smaller than yN. To handle this we will write K = ly/N, and we will let I — > 0. 

Our main purpose in proving this result is to handle i = 2 in (J5.12)) . however, we 
state and prove the result with slightly more generality. For a function g : R i— > R 
define V(x) = V( 9! k,r)(x) to be the function 

V(x) = (g K ( V N (x))-g K (m K (x)))I[m K (x) < R], 

where g K (r] N (x)) = AV\ x „ yl < K/N g N (r)(y)), and g K (p) = E^^^ [g{r] N {y))]\. If the 
configuration rj depends on time, then V will also depend on time, and we shall 
write this as V s . Notice that in the E^ p [g(r] N (y))]} we are considering configurations 
rj on T , and not on the discrete circle. However, this does not change the invariant 
measures, and for this reason we do not alter the notation. 

Lemma 5.3 (Local Equilibrium Principle). Suppose that g is a function satisfying 
\g(x)\ < C(l + x), then 

lim lim sup N 1/2 E NEQ [\ [ < J,V S > ds\] < 0. 

1^0 V^oo|| J || oo<1 J Q 
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Proof. Let X s denote < J,V S >. By the entropy inequality and assumption (HI), we 
have that the following bound for any positive (3, 

+ §. (5.13, 



Eneq[N 1/2 \ £ X,ds\] < ^L\ogE EQ exp{±/3JV j N^X s ds 



It follows from the Feynman-Kac formula that 



E EQ [exp{±(3N [ N^XJs}} < exp{ [ T^ds}, (5.14) 
Jo Jo 

where is the largest eigenvalue of ±PN 3 ^ 2 X S + Ljy. We next reduce the span of 
the Dirichlet form D^^y/f). We have that 

^(V7)>2^i^ f E E E »Mv N (mv x , y x/]} 2 ] 



2K + 



\x-z\<§ \z-y\-- 



where V x>y f = f{j] x,y ) — f{rj). We denote the right hand side of the last line above 
by 2K+1 Djf' x . We may hence bound by 

/V 3 / 2 

iV 3 / 2 sup sup { Eflp [a.V(x).f}-— — rD%°(y/f)}, (5.15) 

/ |a|</3[|J[[oc 2K + I 

because of homogeneity of the system. We next condition on the density of particles 
in (J5.15)) and apply the logarithmic Sobolev inequality of Proposition 14.41 to obtain 
the bound 

/V 3 / 2 

Tf < N 3 / 2 sup sup sup {E UKy [a ■ V(x) ■ f] - C-—H K , y (f)}, (5.16) 

/ \a\</3\\ J\\oo \y\<R A 

for some constant C. Applying the entropy inequality again, we obtain that for any 
M positive the above is bounded by 



(/) (5.17) 



f 1 f i N 3/2 \ 

N 3 ' 2 sup sup sup <^ — log E„ K [exp{aMV}] + — - C-— H K J 

Choosing M such that = C^jp- and letting I = we may combine the last 
bound above with (|5.1H|) to obtain 



ft CtN 1 / 2 C 

E NEQ [N^ 2 \ / X s ds\] < — — sup sup log E UKy [exp{aCl 3 V}] + -. (5.18) 

Jo P l \a\</3\y\<R P 

The next lemma allows us to complete the argument. 

Lemma 5.4 (Local Large Deviations). For any function g which is at most linear 
in k, that is, g(k) < C(l + k), and for all y < R, there exists a constant C = C(R) 
such that for all y 



lo & E i>K,y [ ex P{lAV\ x _ yl < K g(ri(x))}] < C 



K K 2 
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We apply this result to E UK [exp{aCl 3 V}] and obtain that (j5.18|) is bounded by 
C{(3 2 1 2 + r ^/ 2 3 ' 4 } + % fr o m which the result follows. □ 

Proof of Lemma \5.4\ W e begin by fixing a constant Mq > 0. For R > y > Mq/K, 
we proceed as in |CYj . In ()5.19|) . (|5.20|) and (j5.21|) below we look at the single 
site marginal of the grand canonical measures. For simplicity of notation, we omit 
the subscript x from the functions 77. We define g to be the centered version of g, 
g = g-E^\g). 

First, fix 9 so that / (77 - p)e e ^- p ^/ K ^dp p {rf) is equal to zero. A straightforward 
calculation shows that 

9 = -lUM^l + Otf/K 2 ). (5.19) 

K<wn>p 

As we are only considering values of p inside a compact set, these bounds are uniform. 
We use the notation < g; f > p to denote the correlation of g and / under the measure 
p p . Similarly, we define < f\] f 2 ; / 3 > P = Ef Jip [f 1 f 2 f s \. 
Next we define the measure P(rj\9,g,p) by 

dP( V ) = Z-\9, g, p) exp{% - p) + j/ Kg( V )}dp p ( V ), 

with normalizing factor, Z(9,g,p), given by 

J e e ^- p)+ ^ K ~ 9M dp p (r]) = 1 + -i 2 /K 2 H 2 + 0( 7 3 /K 3 ), (5.20) 
where H 2 = < 9; 9 >„. 

Remark 5.5. Note that it is because of this step that we restrict ourselves to linear 
functions g, as otherwise the normalizing function will not be finite. 

We also define the variance under our modified measure: 

a 2 (9,g,p) = Z-\9,g,p)J(r,-p) 2 e^-^ K ^dp Jp (r ] ) 

= a 2 Z- 1 (9,g,p){l-2 1 /KH 1 + 0( 1 2 /K 2 )}, (5.21) 

where a 2 =< 77,77 >„, and i?i = <Th9> j' <r, ' r i^ > f _ <Mh£t£ We now re-write the 
quantity of interest using the notation developed above. 

/-= e e ^\*-y\<K {n ~ p)+1,Km dp p {<n) 
E v K Jexp{lAV lx ^ yl < K g(ri(x))}} = j — — — 

ry-Kia \ P (V = p\0,9,p) ,r QO x 

= Z {9,g,p)--— — — . 5.22 

P(n = p\9 = g = 0,p 
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We now apply Proposition 14.91 with J = 8 to show that 



y/K<T*(6,g,p)P(fj = p\9,g,p) _ 1 + 0{l/R 2y 



^Ka^p)P(fj = p\0 = 9 = O,p) 



We combine this last bound with (j5.19|) through (j5.21j) to show that (|5.22j) behaves 
like 1 + -y/KHi - l/2 7 2 /KH 2 + 0(( 7 + 1 2 )/K 2 ). 

It remains to study the effect of considering g vs. g in the expectation under study. 
Let g* = AV\ x - y \<Kg{j](x)) . By Jensen's inequality we have for any positive a: 

-oT l log E VKy [exv{-ag*}\ < E UKy [g*} < a' 1 \ogE UK y [exp{ag*}}. 

From previous bounds we have 

a' 1 logE VKu [exp{±aAV\ x - y \< K §{T i (x))}] = ±1/KH 1 - l/2a/KH 2 + 0((1 + a)/K 2 ), 

which equals Hi/K + 0(K 2 ) if we select a = 1/K. We combine these results to 
obtain the required result for the case y > 

7-Mog (E UKy [ex P { 1 AV [x ^ yl < K g( V (x))}])= 1 /KH 2 + 0((l + 1 )/K 2 ). 

For y < Mq/K the bound is much simpler as in this case the total number of particles, 
and hence any non-local function, is bounded. In this case the result follows by a 
simple Taylor expansion. □ 

5.2.2. Equivalence of Solutions. By the hydrodynamic scaling limit we know that m K 
and p solve the same limiting differential equations. We use here a standard PDE 
trick to show that the solutions themselves are close. This idea was first developed 
in fHY|. 

Lemma 5.6. Under the assumptions of this chapter we have that 
lim lim E{ [ N- 1/2 V (m K (x) - p(x)) 2 ds] = 0. 

J U xG.T N 

Proof. Let (Af)(x) = f(x + 1) — 2f(x) + f(x — 1) denote the discrete Laplacian, and 
write its inverse as the matrix G = {g(x, y)}. As we work with centered functions on 
the discrete torus, this inverse is well-defined. 

A careful calculation gives that LHY^Hj! is equal to 



L l-^E Y, N ~ 2 9^y)Yt N (x)Y t N (y) 



N 



N < c(r, N ) - (pip), Y»>+J2 <V N (x)) + ^ E <V N (x)). (5.23) 



N 



NONEQUILIBRIUM COLOUR FLUCTUATIONS 



29 



Similarly, we obtain 

L < t - p,Y t N >, 



<r)»-p,c( V N )-<p(p)>. 



(5.24) 



Combining (|5.23|) with (J5.24)) we get 

w^- 1/2 ii^ii;+2<^-p,^>.] 

- E NEQ [N-^\\Y Q N \\l + 2<fj N - p, Y» > 
N- 1/2 H s (x)ds 



E 



NEQ 



Uo 



N~ 1/2 E 



NEQ 



AV x£T nc(t] s (x))ds 



Uo 



(5.25) 



where we let H s (x) = c{rj^{x)) — [c(rj^(x)) — <f(p s (x))}{r]^ [ (x) — p s (x)} 
We next replace H(x) with its local average 



H K (x) = AVy 



\x-y\<KjN 



c(v(y)) - {c(v(v)) - <p(p(x))}{v(y) - p{x)} 



Assumptions (LG) and (M) imply that we have finite exponential moments for the 
measure p p . By the entropy inequality and assumption (HI) we conclude that 



N- 1 ' 2 E NEQ [ f AV x ^c{^{x))ds] = OiN- 1 ' 2 ). 
Jo 

Thus, using summation by parts, along with I = we obtain that 

Jo 



(5.26) 



NEQ 



0(1). 



Combining the above result with 1)5.25)1 we have 

Eneq [N~ ^ 2 1 1 Y t N 1 1 1 +2<fj-p,Y t N >„] - E[N-W\\Y»\\l + 2<fj-p, Y N >„] 



E 



NEQ 



f N^ 2 Y,H?{x)ds 

J ° xdf N 



+ 0{1) + 0{N 



-X/2\ 



(5.27) 



We will make use of the local equilibrium principle to handle the remaining term. 
Define F(r)(y), p(x)) to be the quantity 

c(v(y)) - {c(v(y)) - <p(p(x))}{n(y) - p( x )}, 

and let F(m K (x), p(x)) denote its expectation under the measure p p=m K {x)- A brief 
calculation shows that F is equal to I ip{m K (x)) — <p(p(x)) \ \ m K (x) — p{x)) \ . Lastly, 
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let s K (x) = AV| a ._j / | < k rf(x). We split H K (x) into three pieces: 



H K (x) = AY, 



\x-y\< 



F(rj(y),p(x)) 

-F(m K {x),p{x))] I(m K {x) < RDs K (x) < R 2 ) 
+ AV ]x _ yl < § F(r)(y),p(x)) l(m K (x) >RUs K (x) > R 2 ) 
+ AV {x _ yl <K F(m K (x),p(x)) l{m K {x) <Rns K (x) < R 2 ) 



AY 



\x-y\< 



fW(x)+fW(x)+F®( 



x) 



By the local equilibrium principle the term involving F^ vanishes in the limit. Even 
though the term F — F is not linear in r] as required in the local equilibrium prin- 
ciple, the addition of the condition I(s K (x) < R 2 ) fixes the problem pointed out in 
Remark 15.51 . Applying bounds from Remark 13.21 of Section l3~2l we have that 

AV ly _ xl <KF(ri(y),p(x)) < c m K (x) - c.AV^K^y) - p(x))\ (5.28) 

for some constant c > 0. If m K (x) > R, we fix C and choose R so that 2| |yo||oo < R—C, 
which implies that m K (x) — p(x) > C and hence 1)5.28)1 is smaller than 

—A (m(x) — p{x)) 2 I(m K \x) > R) 

We next pick C large enough so that A = c\ — 2c/ C > 0. 

If, on the other hand, m K (x) < R, then we control (j5.28j) using s (x) > R 2 - For 
R 2 large enough 

{r 1 {y)-p{x)f>R 2 -2R\\ 

P\\oo- 

Hence, for c > take R 2 > R/c + 2i?||p|| 0O , we get that R < c(R 2 — 2||p|| 0O i2) 
implying 

cm K (x) <cR< ccAV ly _ xl < K/N (r](y) - p(x)f . 
We thus obtain a bound on (|5.28|) of 

-B(m K (x) - p(x)) 2 I(m K (x) < R,mg(x) > R 2 ) 

where B = c\ — cc is positive as long as c was chosen to be sufficiently small. Thus, 
we obtain that for fixed R and R 2 large enough, and A = min(A, B) > 

AV\ x _ y \< K/N F {2 \x) < ~A(m K (x) - p(x)fl(m K \x) < R, mf (x) > R 2 ) (5.29) 

A similar bound holds for F^\ using the fact that <p' > C\. Combining these last 
bounds together with ()5.27|) we find 

E NEQ [N-^\\Y t N \\l + 2<fj-p, Y t N >„] - E[N-^\\Y N \\l + 2<fj-p, Y N >J 
< -A*E NE q[ [ N~ 1/2 (m K {x) - p{x)) 2 l{m K \x) > R U mf (x) > R 2 )ds] + o(l) 
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where A* = A + c\. We shall now use assumption (Fl) together with the inequalities 

1 

2' 



±2 <fj-p,Y t N >*< 2N 1,2 \\f]-p\\l + ]-N- 1,2 \\Y t N \\l (5.30) 



and | \r] — p\ \ 2 < C\\r) — p\\ 2 , to conclude that 

\N-^E NEQ [\\Y t N \\^ < -cE NEQ [ f N- 1 ' 2 J2 (m K (x)-p(x)) 2 ds}+o(l). 

The result follows. □ 

5.2.3. Conclusion. By Lemmas 15 .31 and 15. 61 we have that (j5.12J) holds for i = 2,3. We 
consider Q± next. By summation by parts and (j5.26|) we have for / e C 1 (T) 

\E[f <f,*?(s)>ds]\<Cl, (5.31) 
Jo 

Using arguments similar to those used to obtain ()5.29|) we show that for i — 4,5 

| < 1, &?(x) > | < C'{m K (x) - p(x)) 2 l[mf (x) > R] (5.32) 

for some positive constant C and then apply Lemma [5.61 This concludes the proof 
of Theorem 15.21 

5.3. Quadratic Variation. We need to show that Y t has quadratic variation given 
by (|5.4|) . We know that under P^, N^(f) defined in (|5.6J) is a martingale. Hence, it 
remains to show the following proposition. Note that it implies that test function / 
must be at least C 2 (T). 

Proposition 5.7. For a function h in C 1 (T). 



E 



< h , c(Vs) ~ <P(Ps) > ds 



0. (5.33) 



Proof. This follows from the Boltzmann-Gibbs principle of Theorem 15.21 which is 
a stronger result, in combination with the hydrodynamic scaling limit of Proposi- 
tion EM □ 

5.4. Unique Solution to the Martingale Problem. We recall here a result due 
to Holley and Stroock |HSj which guarantees the existence of a weakly unique solution 
to our martingale problem. 

Let B t denote the nonpositive operator |D(p t )A and V D {t, s) the associated semi- 
group. Also, let St be the operator a/ (p(p t )V. 

Theorem 5.8. Fix a positive integer m > 2. Let P be a probability measure on the 
space {(C[0, T], H_ m ), J 7 }, where T = U t >oJF t; and T t is the canonical filtration of 
the process <Y t ,f>. Assume that for any f G C°°(T) 

<M t ,f> = <Y t J>-<Y ,f>- f*<Y 8 ,B a f>ds . . 

and <M t ,f> 2 -f*\\S s f\\ 2 ds ^ ) 
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are L l (P) martingales with respect to Tt- Then, for all < s < t, f in C°°(T), and 
subsets A ofM, 



P[<Y t J>eA\F a 



(y-<Y s ,V D (t,s)f>) K 



f s \\S u V D {t,u)f\\ldu I 2 L \\S u V D (t,u)f\\ 2 2 du 

P almost surely. In particular, \5. 34)) together with a unique initial distribution for 
< Y , f > uniquely determines P on {(C[0, T], H- m ), J 7 }. 

Thus, to guarantee uniqueness of the limiting measure P, we need only check that 
we have a unique initial measure; this is simply condition (F2). 

5.5. Tightness. Let Pjy denote the probability measure on D([0,T],TC- m ) induced 
by the fluctuation field Y t N started in the measure (q , and E NE q the associated 
expectation. In this section we prove that the sequence of probability measures P^ 
is tight. A consequence of the proof will be that the limit points are concentrated on 
the space of continuous paths C([0, T], TC- m ). 

For a function F in D([0, T], W_ m ), define the (uniform) modulus of continuity, 
ujs{F), for a fixed 5 > 0: 

u s (F) = sup \\F t -F s \\_ m . 

\s-t\<S,0<s,t<T 

To simplify notation we will often denote this supremum simply as sup s t \ \Ff — F 8 \\- m . 
By the well-known Arzela-Ascoli result, it follows that to show the sequence Pjv is 
tight we need to show that it satisfies two conditions: 

(Tl). lim M ->oolinisup Jv ^ oo Pjv[sup < t < r ||y t ||_ m > M] = 
(T2). lim^o limsupjv^ P N [u s {Y) > e] = Ve > 0. 

We will also make use of the result, due to Aldous, |Aldj . For a proof, see for 
example Proposition 1.6, of Section 4 in |KL| . We will state it in some generality. Let 
X denote a Polish space with metric d(-, •). For a function g in D([0,T], X), define 
the modified modulus of continuity, uj' s (g), for a fixed 5 > 0: 

^a(#)=inf max sup d(g t ,g s ), 

{ti}V =0 0<Kr ti < s<t< t i+1 

where the infimum is taken over all partitions {^}[ = o of [0, T] such 
that t = < ti < . . . < t r = T and t i+ i —ti>8. 

Proposition 5.9. A sequence of probability measures Pjy on D([0,T], X) satisfies the 
condition 

limlimsup Pn[uj' s (Y) > e] = 
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provided that 



lim limsup sup P^[d(Y T , Y T+e ) > e] = 0, 



(5.35) 



where T denotes the collection of all stopping times bounded by T . 



Proposition 5.10. The sequence of measures Pn is tight in D([0,T],H- m ) , for any 
m > 4. Moreover, all limit points are concentrated on continuous paths. 



We first consider condition (T2). Recall that Y S N is the field vN(c(rj^ (x) — p t {x)). 
We may hence write 



rN 71 <fN I 



\Y t N - Y s N \\. m < sup | < A N f, / Y u N du > | + ||Mf - Mi „_,„ 



where < M f , / > is a martingale with quadratic variation given by the martingale 
in ()5.6|) . Using this result, we reduce the proof of condition (T2) to showing that 
both quantities 



(T2), 



Pn sup^supu/u^! | < A N f,f*Y*du > | > 
(T2) B Pn [sup S) , ||Mf - Mf ||_ m > §]. 
decrease as 1 1 — s | < <5 converges to zero. We begin with (T2) B , and we split the proof 
into several lemmas. 

Lemma 5.11 (Bounding the tails 1.). There exists a finite constant C such that for 
every eigenfunction e z ofH m , 



lim sup E NEQ 

N^oo 



sup | < Mf , e z > 

0<t<T 



< CT{1+ < Ae z ,Ae z >} 
Lemma 5.12 (Bounding the tails 2.). For m > 3, 



lim lim sup E NE q 



sup 

0<t<T 



J2i<M t N ,e 2 >}\ 



2 -m 



;>K 



0. 



This last lemma implies that to prove condition (T2)b is satisfied, we need only 
check that for bounded K 



lim lim sup Pn 

N^oo 



sup i< M t N - M?, *. >} 2 7," m > ^ 



\z\<K 



0. 



This will follow if we can prove the following result. 
Lemma 5.13. For any f in C 3 (T) 



lim lim sup P/v 

<5^0 at^oo 



sup|<Mf -Mf,/>|>e' 

s.t 



0. 
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We make use of the following fact 

u s (< M t N ,f>) < 2cu' s (< M t N ,f >) +sup| < M t N ,f> - < M t N J> |. 



By the definition of M t we know that 



| < M?,f> - < M{1J> | = | < Y»,f > - < Y«,f> |, 

and notice that because at most one particle makes a jump at any given time, this 
quantity is bounded above by 1 |oo-^ -3 ^ 2 - We therefore obtain that 

lus(< M t N J>) < 2u' s {< M t N J>) + HfllooiV- 3 / 2 (5.36) 

Proof of Lemma \5.1S[ By ()5.3fij) above and Proposition 15. 91 it is enough to check that 

lim lim sup sup P N [\ < M? +e - , / > | > e'l = 0, (5.37) 

e<8 

By Chebychev's inequality, and because both r and 9 are bounded stopping times, 
we may bound P N [\ < M? +9 - Mf , / > | > e'] by 



^E NEQ 



l<V?(x)WNf(x) 2 + V N f(x + l/Nf]ds 



□ 



Applying the results of Section 15.31 completes the proof. 

We finish the argument for (T2)b by proving Lemmas 15.111 and 15.121 

Proof of Lemma \5.11\ Because M t is a martingale we may use Doob's inequality to 
bound E NE q[sup 0<t<T \ < Mf ,e z > | 2 ] by 4E NE q[\ < M^,e z > \ 2 ] and this is equal 
to 



4 
2 



[ [c(^(x))][V N e z (x) 2 + V N e z (x - 1/N) 2 + V N e z (x + l/N) 2 ]ds. 

Using the convergence established in section 15. 3^ and the fact that <p is bounded 
uniformly, we may bound the above by T||<^|| 00 {1+ < A N e z , A N e z >}. □ 

Proof of Lemma \5.1 6 A This follows since by Lemma ()5.11|) the quantity 



E 



NEQ 



sup 

0<t<T 



£{<Mf,e 2 >} 2 7 r m 



N„ A N, 



is bounded by C£ N >x7r m {l+ < A N e z ,A N e z >}. 
We next turn to condition (T2)a- 



□ 



(T2)a lim lim sup P/v 

<5^0 7V-+oo 



sup sup 

*J ll/llm<l 



<Y u N ,A N f>du 



e 

>2 
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where Yj* is the field \/~N {c{j]f{x) — <p(pt(x))). Similarly to the proof of the Boltzmann 
Gibbs principle, we split this field up as the sum of the following terms: 

K = VN{c K {r ] f{x))- V {mf{x))}l[mf{x)<R] 
tf* = VNif l (p t (x)){mf(x)- Pt (x)}I[mf(x)<R] 

VN{ip{mf(x)) - <p(p t (x)) - <p'{p t {x)){m?{x) - p t (x)}}l[m? (x) < R 
-VNcp{ Pt {x))I[m?{x) > R] 
ip*f t = y/Nc K (x)I[m?{x) > R] 



3,t 



5.38) 



It hence remains to prove that for i 



lim lim sup P/v 

5^0 iV-*oo 



sup sup 

M ll/llm<l 



1,...,6 

t 



<*»A N f>du 



> 



12 







For i 2, 3 we first apply Chebychev's inequality and use Lemma [5.91 For i — 1 we 
argue as in (j5.31|) . For i = 4 we have Lemma f5. 61 Lastly, we use (j5.32j) and Lemma 
15.61 to handle i — 5, 6. 



P, 



N 



sup 



1 <1 



T+6 



<*»A N f>du 



> e 



< -E 



NEQ 



sup 

f\\m<U0 



T 



<*»A N f>\du 



It remains to handle the terms with i = 2,3. In order to do this we use a special case 
of the well-known Garcia-Rodemich-Rumsey inequality. 

Lemma 5.14 (Garcia-Rodemich-Rumsey inequality). Given a function g and a strictly 
increasing function ip such that ^(0) = and lim^oo ip{u) = oo, then 



\g(t)-g(s)\< 



We choose ip(x) 
obtain 



where 

_ e N\x\ 



B 



o jo 



\9{t)-g(s)\ 



u 



)dtds. 



1 in the above inequality. By integration by parts we 



M i+ ^)^* s ^ {iog(i+ ^ )+4} ' 

To simplify notation we drop the subscripts on the field tyf s and simply write it as 
^/(s). We apply the Garcia-Rodemich-Rumsey inequality to any such field ^(s) in 
the following way: choose g(t) = J* < ^(u),f > du, where ||/|| m < 1. We then 
obtain 

| J < (u), / > du\ < V~5{\og{\ + ^) + 4}, 



.'!(:» 



HANNA K. JANKOWSKI 



and we may also bound B by 



'o jo 

and we have the inequality 



exp | Ar M ^WH U^ J _ 1 j > ^ 



A = sup|| / ^( M )^||_ m <i" 1/2 {log(l + 45 r 2 )+4}, 



(5.39) 



where 



B = J J exp I AT s j dsdt - T . (5.40) 

It follows that 

5o< je^ 16 ^ ^ < e /2 + ^, 

and this last quantity is smaller than e for large enough N. We thus have that for 
any field \I> the following inequality holds 

P N (A >e)<P N (b > j e»<fl*^ . (5.41) 

Next, we apply the entropy inequality (see for example Prop. 8.2, of Appendix A 
in |KLp to bring the problem back into equilibrium measure. We also use here 
assumption (HI). We thus obtain that for any set A such that P^q(A) < e - NC ( s ) 

P N [A] < [log 2 + ^(0)][log(l + PgaiA)- 1 )]- 1 , 

< [CN)[\og(l + P» Q (A)-i))-\ 

< l/C{5) (5.42) 

where Peq is the measure induced by the field Y N started in the equilibrium measure. 
We have thus reduced the problem to showing 



f Bo >^e N ^A <e' NC ^, (5.43) 
where C(5) — > oo as 5 — > 0. We now apply Chebychev's inequality 

P&, (ft > ^exp {N±f*}) < 4»lexp {-N±S^} . 
Suppose next that we have that 

Eeq[B q ] <e CN , (5.44) 
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where C is some positive constant. Notice that x < e~ x for small x, and hence if 
(|5.44j) holds, we have that there exists a <5 > and a positive constant C depending 
on 5q such that for all 5 < 5q we have 



Pg, (Bo > £ exp {N^}) < exp [-NC^ _ ^r^} . 
Hence (I5.43J1 holds. We have in fact proved the following. 



(5.45) 



Proposition 5.15. Suppose that there exists a constant C > and an N such that 
for all N > N 



E 



EQ 



JO 



exp N- 



f s *(u)du\ 



Vt 



dsdt] < e 



CN 



(5.46) 



then 



lim lim sup Pm 

6^0 N 



sup 1 1 / ty(u)du\ |_ m > e 



0. 



Hence it remains to prove (|5.46|) . We follow a similar argument presented in |CYj . 
By the Cauchy-Schwarz inequality 

ft n 2 N T l N 

N\\ / m u du\\_ m {t-s)- l / 2 < —|| / ^ w || 2 _ m + - 

t s J s a 

and hence to obtain ()5.46j) it suffices to prove that for a constant C = C(a) > 0, 
independent of i and s, 

-t 



E 



EQ 



exp < N 



(t-s) 



®udu\\ 



< e 



CN 



(5.47) 



Define Q to be a Gaussian measure defined on 7i a with zero mean and covari- 
ance S = (1 - A Ar ) _/3 , where A w / = iV 2 [/(a; + 1/iV) - 2f(x) + /(x - 1/N)}. The 
existence of such a measure is guaranteed for (3 > 1/2 by the theory of Gauss- 
ian measures on Hilbert spaces (see for example jVarj). Define the inner product 
< F,G > k =< F, (1 - A N ) k G >. We thus have that 



exp{< F, G >}dQ(G) 



exp{< F, (1 - A N y a G > a }dQ(G) 
1 



exp 



-a-/3 



We choose m = a + (3 and changing the order of integration we may write (|5.47j) as 



E 



EQ 



exp 



aN l l 2 (t 



-1/2 



< / $(u)du, f > 



dQ{f) < e 



C{a)N 



(5.46 
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Suppose that we could show the following bound for some positive integer k 



exp jaiV 1/2 (t- S )- 1/2 < tf(u)dit,/>J 



< e Ca 2 {||/|| fe +l}_ (5 _ 4g) 



Plugging condition ()5.49|) into (|5.48j) we would obtain that the left hand side of (|5.48jl 
is bounded by 

exp{Ca 2 } J exp{Ca 2 <f,f > k }dQ(f), (5.50) 

which, because Q is a Gaussian measure, integrates to 

{det (1 - 2Ca 2 (l - A N ) k - m )}- 1/2 , 

for a sufficiently small and m — k > ~. Since this is (quite) smaller than e CN for some 
positive constant C we have reduced the proof of ()5.47|) again to showing (|5.49jl . That 
is, we have proved the following result. 

Proposition 5.16. Suppose that there exists a constant C > and an No such that 
for all N > N 

E EQ [exp{aN 1/2 {t - s)~ 1/2 < f, ^(u)du >}] < e Ca2{m " +1} . (5.51) 
then for m — k > | 

E EQ l£ £ex P dBdt] < e c». 

We thus need to show that (|5.51|) holds for ^f t when i = 2, 3. To do this we 
will repeat the arguments of the local equilibrium principle for both terms. Again, 
to simplify notation, we will denote simply as Both terms satisfy the same 
bounds. As before, we begin with the Feynman-Kac inequality, 

\ogE EQ [e^<^n^>] < SN^snp sup (^JtII/IU^)^] - N^D K>y (^f)\ . 

h \y\<R { K ) 

We apply the logarithmic Sobolev inequality given in Theorem 14.41 and the entropy 
inequality as in (j5.16|) and (j5.18|) to bound the line above by 

f /V 5 / 2 1 
6N ^ sup i og £ [exp{||/|| oo7 ir 3 /iV 5 / 2 v[/(a;)}] , 

where 7 = ^_ s - We next apply Lemma 15.41 and obtain that this is smaller than 

C(l, Nq, cio)(||/||oo+ I \f\ lie) * 2 f° r a ^ a > a 0) an d N > N . Because we are on the unit 
circle T we have, using the Sobolev inequality sup x&T f(x) < J T \f'(x)\dx, for some 
new constant C{1) 

XogEEQ^iaN^it-sY 1 / 2 <A N f, f *(u)du>}\ < C(l)a 2 {l + \ \f\ | 3 }. 

J s 
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Combining this with Remark 15 . 1 51 and Remark 15.161 we may conclude that for % = 2, 3 



E EQ [J J exp(iV(t 



-1/2 I 



^i(u)du\\- m J dsdt] 



< e 



CN 



(5.52) 



This concludes the proof of condition (T2). 

By Chebychev's inequality as well assumption (F2), to prove that condition (Tl) 
is satisfied it is enough to show that 



lim 



M-*oo Pn 



sup t sup|| / || m < 1 {< A N f,f*Y^du >}> M 
and E NE q [sup t || || _ m ] < oo 
We begin with the latter. By Lemma 15.121 it is enough to prove that 







(5.53) 
(5.54) 



E 



NEQ 



sup<Mf,/> 



< oo. 



This follows from Doob's inequality and the result of Section I5~3l 

The case (I5.53J1 will be a repeat of the argument for (T2)a- We split the field Y as 
in (|5.38|) . and show that 



lim 



M->oo r N 



Pa 



sup t sup | | /||m < 1 {< A N f, f* *»du >}> M 



0. 



for % = 1, . . . , 6. For i ^ 2, 3 we first apply Chebychev's inequality. For i = 1 we argue 
as in (|5.3ip . For i = 4 we use Lemma l5~o1 Lastly, we use (|5.32|) and Lemma I5~o1 to 
handle i — 5,6. For the two remaining cases, it turns out that we have already done 
all of the work. We again use the Garcia-Rodemich-Rumsey inequality, Proposition 
15.141 We need simply to change the roles of S and e. Again, we first argue for a 
general field Let A and B be as in (J5~33|) and (l5~47IJ) . From (Jo~4"T| and (l5~l2l it 
follows that for any field if we can show that 

T 2 



Ne/WVT^ < 



-NC(e) 



where C(e) 



oo as e — >• oo, then 

lim^oo lim sup^ P^r [sup t 1 1 J Q * *ff(u)du\ 



> e 



0. 



(5.55) 

By Chebychev's inequality, it remains to show that Eeq[Bq] < e CN for some positive 
and finite constant C. However, this is exactly the content of (|5.52|) . This concludes 
the proof of tightness. Notice again that we have in fact proved the following. 

Proposition 5.17. Suppose that Eeq[B ] < e CN , with B defined in l[5.40i) , for some 
positive and finite constant C , then 

-t 



lim lim sup Pjv 



sup 

t 



^{u)du\\- m > e 



0. 
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6. Colour Density Fluctuations 

In this section we prove Theorem 11.11 For the sake of brevity, we shall write out 
the proof for the case when k = 2; the generalization to any value of k is immediate. 
As for the colour-blind key ingredient in the proof is the logarithmic Sobolev 

inequality. We use the additional condition (E) in order to have the logarithmic 
Sobolev inequality, and the remainder of the work is valid without this assumption. 

The general outline of the proof is the same as the approach of the previous section. 
We show that the sequence P^ is tight. We then analyze the martingales under 
possible weak limits of Pn and use them to identify the form of the martingales 
under the limiting measure. We finish by showing that there is only one measure 
which solves the martingale problem. 

It is simpler to work with the following version of ()1.6|) . for any smooth test func- 
tions fi 

< dY iit , U > = < Y iit , djipif > dt+ < Y jit , dpitpj" > dt+ < dWl yfoft > (6.1) 

We begin by identifying the martingales. For / in iff define < Y t N , f > to the 
pair {<Y t N 1 ,f 1 >,< Y t N 2 , f 2 >}- Under P N we have that 

< M t N ,f>=< Y t N ,f> - < Y N ,f>-f* < Y S N , Af > ds (6.2) 

is a martingale, where Y denotes the coupled fluctuation field created by the pair of 
fields VN ( Ci (r/(-)) - <pM-))). 

The outline of this section is as follows. In Section I6~T1 we identify the drift martin- 
gale under the limit of any converging subsequence of the measures Pn- To do this 
we prove the Boltzmann-Gibbs principle. We next handle the quadratic variation and 
show that there is a unique solution to the martingale problem. We finish the section 
by proving that the sequence of measures Pn is tight. As in the previous section this 
completes the proof of Theorem 11.11 

6.1. Identifying the Drift. Arguing as in the previous section in (|5.7jl and (|5.8jl . 

to show that the limiting martingales ()6.2j) are the same as the drift martingale of 
(|6.1|) . it is enough to prove the following. 

Theorem 6.1 (Boltzmann-Gibbs). For any f in C 1 (T) and i = 1,2 



lim sup E NEQ 

N—foo 



[ < f, VN{ci{ri?) - ipi{p B ) - V^(p s ) • (r?f - Ps)) > ds 
Jo 



0. 



The first step of the proof is to replace Pt(x), the weak solution of d t p = |V '£> 2 (p) V p 
with the solution to the discretized version 

d t p = ~V N D 2 {p) VjvP (6.3) 

with initial conditions Pt(x)\t=o — Po{ x ) f° r x m ^ N x ^ N • The difference between the 
two solutions is of order [RM. and hence does not affect our calculations. Because 
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of this fact and in order to simplify notation, we continue to denote the solution of 
(|6.3|) as p in the remainder of this section. 

We proceed with a proof of the two main ingredients: the local equilibrium principle 
and "equivalence of solutions". Notice the similarity to the proofs of the previous 
section. This happens because the jump rates as well as their expectations for the 
individual types are of the form <pi = p l ^p, exhibiting linear behaviour only in the 
density p\ 

6.1.1. Local Equilibrium Principle. For a function g : N + x N + i— > IR define 

V(x) = V {g ^ R) {x) = AV ly . xl < K/N {g(v N (y)) - E^ mK Jg(r, N (y))]}l[m K (x) < R] 

= {9K{ri N {x)) - g K {m K {x))}l[m K {x) < R], 
with m K \x) defined as in the previous section. 

Lemma 6.2. Suppose that g is a function satisfying \g(ki,k 2 )\ < C(l + ki), for all 
(ki, k 2 ) and for some i, then 

lim lim sup N 1/2 E NEQ [\ [ < J,V S >ds\] < 0. 

i^0 7V^oo|| J || oo<1 J Q 

Without loss of generality, we will consider only the case g(ki, k 2 ) < C{1 + k\) in the 
proof. We begin with a lemma. As in the previous section we state the next result 
for configurations 1] prior to the diffusive re-scaling. For g as in Lemma I6.21 define 

9 = 9- E ^i =m K (x)ti=lt2 {-\v 2 )[9}, where mf(x) = AV\ x _ y \< K ^{x) . 

Lemma 6.3 (Inhomogeneous Local Large Deviations). There exists a positive con- 
stant C independent of the configuration i] 2 such that 

logE VKy( . M [exp {jAV^y^Kgirtix))}] < C{^ + 



Proof. We first note that r\\ conditioned on r] 2 satisfies the conditions of a nonho- 
mogeneous zero-range model as was shown in Section 14.11 and hence we have the 
necessary local limit result from Theorem 14.101 with bounds which are independent 
of 7/2 (and hence the site x). We then repeat the same argument as in the proof of 
Lemma 15.41 □ 

Proof of Lemma \6. 6 A Let X(s) denote < J, V s >, as before. Most of the proof is the 
same here as in the colour-blind case, and hence we omit many of the details. By the 
entropy inequality and Feynman-Kac formula, we have 

E NEQ [N 1/2 \ f X{s)ds\] < -^-tsnpT N (s) + (6.4) 

JO P ly s P 
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where T N (s) is the largest eigenvalue of ±PN 3 ^ 2 X(s) + L^. The biggest difference in 
this setting is that we need to introduce an additional conditioning. We reduce the 
span of the Dirichlet form as in ()5.15|) and condition on the average density and the 
particles of colour two to obtain that T N (s) is bounded above by 

/ |a|<^||J||ac \y\<R m zi\ -t- i 

where -D^ (/) defined in ()4.3|) . Using the logarithmic Sobolev inequality (The- 
orem 14. 5 j) , the entropy inequality and Lemma (jfi.Hj) as in the colour-blind proofs we 
conclude that 

E NEQ [N^\ j\(s)ds\] < C{(3 2 1 2 + + |, 

from which the result follows. □ 

6.1.2. Equivalence of Solutions. Define to denote the centered version of the func- 
tion Y/J as in (|3.1U|) . Because of the symmetry of the process, we only consider | \Y^ t \ | 2 . 
The same formula holds as in the colour-blind case: 

EneqIN-^WY^II +2<fj»- h , ?» >.] 

- E NEQ [N-W\\Y»\\l + 2 < fj» - p u Y» >.] 

= E NEQ [ f{- < Cl (r, N ) - Vl {p)X N s > +iV~ 1/2 (l + 1/N) c^mds] 



E NEQ [ [ N^AV^NH^ds] + N-V 2 E[ [ AV x &»c x {rfi{x))ds] 
Jo Jo 



(6.5) 



where H\(x) = ci(r] N (x)) — (ci(r) N (x)) — ipi(p( x )))(Vi ( x ) ~ P X ( X ))- 
We next replace H\(x) with 

H?(x) = AV\ x -y\< K /N[ci(r] N (y)) - {dfa^y)) - V i(p(x))}{^ (y) ~ A*)}]- 

By integration by parts and (J5.26)) we conclude that the error produced by this re- 
placement is of the order of 0(1). Line (|6.5|) is thus equal to 

Eneq[ [ N^ 2 AV xeT NH^ s (x)ds] + 0(1) + 0(N- 1 ' 2 ). (6.6) 
Jo 

Next define F 1 (r l N (y), p(x)) = Cl (ri N (y)) - {c^ (y)) - Vl (p(x))}{^ (y) - p^x)}, 
and split H K (x) into three pieces: 

Rf = AV lx . yl < K/N [E 1 (r } N (y),p(x))-E 1 (m K (x),p(x))]I(m K < R)I(sf (x) < R,) 
+AV lx - yl < K/N [F 1 (r t N (y), p(x))]I(m K > R U sf (x) > R x ) 
+AV lx „ yl < K/N {F 1 (m K (x),p(x))]I(m K < R)I(sf(x) < R,) 
= AV\ x . y \< K/N [F[ l) (x) + F?\x) + FP(x)}, 
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where 

F 1 (m K (x),p(x)) = E„ p=mK jF 1 (r 1 N (y),p(x))} 

= -(Mm K (x)) - vi(p(x)))(m?(x) - p\x))], 

and sf(x) = AV\ x - y \< K / N {r/f (x)} 2 . 

By the local equilibrium principle the first term vanishes in the limit. Even though 
the term F is not linear in t] as required in the local equilibrium principle, the addition 
of the condition I(sf(x) < Ri) fixes the problem originally pointed out in Remark 
15.51 By an identical argument to the one presented in the previous section, for fixed 
R and R\ large enough there exists a constant A 2 > such that the following bound 
holds: 

AV\ x - v \< K/N Fl 2 \x) < -A 2 (m? (x) - p 1 {x)fl{m K \x) > R U sf (x) > Rj. (6.7) 

We handle F{ similarly. In this case we make use of the fact that the function ip\ (p) 
is strictly increasing (Proposition 14.3)1 . We thus have that there exists an A3 > such 
that: 

AV ]x ^< K/N F^{x) < -I 3 (mf (x) - p\x)) 2 l{m K {x) < R)I(sf(x) < R x ). (6.8) 
Combining the arguments from (J6.6j) through (|6.8|) . we find that 
E NEQ [N~ 1 / 2 \\Y t N l \\l + 2 < ft - - Pl X N i >*] " ExeqIN-^WY^WI + 2<f} 1 - p^Y* >.] 
< -A 2 A A 3 E NEQ [ f N- 1 ' 2 (™?(x) - p\x)fds] + 0(1-') + o(l) (6.9) 

Using the bound (|5.3U|) together with assumptions (Fl) on ()6.9|) as in the colour-blind 
proof, we conclude 

Theorem 6.4. 

lim lim E NEQ [ [ N- 1/2 V (mf(x)-p i (x)) 2 ds] =0, i = 1,2. 

We are now ready to prove Theorem 16.11 Define the field <& N (t) = {3>x (t), <3>^ (t)}, 
with $>i(t) given by 

{ci(T/f ) - cpxipt) - dxcpxipt)^ - p 1>t ) - d 2 (pi{pt)(r]2 7t - p 2 ,t)} , 

and Q 2 {t) defined similarly. 

As before, we split <& N into five separate parts. 

& N (t) = *»(t) + *?(t) + &%(t) + <frf (t) + (t), 
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where the first entry in each is given by 

IT 1 ^ = Ci (t^(*)) -cf (r^(:r)) 

-d 1 ip 1 {p{x)){r]^ {x) - mf (x)} - d 2 <pi(p(x){r)% (x) - mf (x)} 

N- 1/2 ^ 2 = {cf (r} N (x)) - Vl {m K {x))}l[m K {x) < R] 

N- 1 "^ = [<pi(m K (x)) - ipi{p{x)) 

-9iV9i(p(x)){mf (x) - p\x)} 
-d2<pi(p(x)){m$(x) - p 2 (x)}] I[m K (x) < R] 

N- l 'V 1A = \-^{p{x)) 

-divx(p(x)){m%(x) - p 1 {x))} 
-d 2 Vi{p{x)){m%(x) - p 2 (x)}] I[m K (x) > R] 

ATV2$iv = E ?(x)I[m K (x) > R], 

with cf-(x) = AV\ ^KCiirjijj)) . By symmetry, we only work with the &i(t). By 
the local equilibrium principle and Theorem 16.41 we already have for i = 2 and 3 that 

EneqII I <$&,/><fc|]-0, (6.10) 
Jo 

for any continuous function /. The result for i = 5 follows from the same argument 
for the colour blind case, because cf(x) < c K (x). Next, we consider i = 4. We may 
pick R sufficiently large so that 

N 

\^(t)\<CN^ 2 Y,(™«(x)-p l (x)) 2 

1=1,2 

and hence (J6.10)) holds also for i = 4. It remains to consider i — 1. An argument using 
summation by parts like the one presented in the previous section proves the result. 
Notice that it is in this part of the argument that we require the most smoothness of 
the function /. It is enough to consider test functions / which have a bounded first 
derivative. This concludes the proof of the Boltzmann-Gibbs principle. 

6.2. Quadratic Variation. We begin with a lemma. Let K « N, fix R large, and 
define 

X*(x) = AV ]x _ y] ^{ Ci (^(x)) - Vi {p t {x))}l[m K {x) < R]. 

Lemma 6.5. Assume that \t — s\ < T. For any choice of I = ^= there exists a 
constant C = C(l,R) so that the following bound holds uniformly for ||J||oc < 1- 

log£ BQ [exp{ 7 iV < J, / X*du >}} < C( 7 2 + l)T. 
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The same bound holds if we replace J by a function dependent on time with uniform 
bound \\Jt\\oc < 1- 

Proof. By the Feynman-Kac inequality, and conditioning as before, we have that the 
left hand side is bounded above by 

TN sup {E^ [yXifl] - ND NtP (y/f)} 

/>0:£ Vip2 [/]=1 

< TN sup sup sup sap{E VKtvi .\ rKt) [yXi{x)f]] - ^^^(.^(w)} 

f x \y\<R r? 2 71 

We next apply the logarithmic Sobolev inequality, Theorem 14.51 followed by the 
entropy inequality to obtain that this is less than 

iV 3 

T— sup sup{log^ K!/( .|, )2) [exp{7K 3 /iV 2 X i }]}. 
A \y\<R m 

Using Lemma in the above along with ly/N = K proves the result. □ 

We now turn to the calculation of the quadratic variation. For fixed ct\ and we 
have that a\ < M£[, fi > +c*2 < M^ 2 , /2 > has quadratic variation under P/v given 
by: 

Ey f N ~" Y,^(x)) + c t (v?(x + ^W N f(x)} 2 ds. (6.11) 

To show that in the limit these martingales match with (jfj.ljl it is enough to prove 
the following proposition holds for i — 1, 2. 

Proposition 6.6. For a smooth function h in C 1 (T). 
Eneq 

Proof. Let X Sji denote the field c^rj^) — ipi(p s ) for i — 1,2. Our first step is to replace 
X S)i with an average over boxes of size K. By summation by parts, the replacement 
adds a term of size |r. Hence we need only worry about the inequality using the field 
averaged over boxes of size K : Xf t . Fix a large R and consider 

X*(x) = AV\ x . y \< K/N [airfHx)) - wMx))] l[m K (x) < R] 
+^V\x- y \<K/N [ci(n?(x)) - (pi(fi a (x))] l[m K {x) > R] 
= X^\x) + xff\x) 
By Chebychev's inequality the second field, Xf'^ b \x), is bounded by 

ENE Q [j\<KX^ b) >\ds]<^ 



< h } d(ri s ) - (fi(ps) > ds 
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for some constant C. By the entropy inequality, it follows from Lemma Id 51 and as- 
sumption (HI) that 

cp C 



E NEQ [[ <h,X^ a) >ds}< 
Jo 



N (3 

The result follows if we first let iV converge to oo, then do the same for /3, then R. □ 

6.3. Uniqueness of the Martingale Problem. In this section we show that the 
martingale problem described by the generalized stochastic differential equation (jl.6|) 
has a unique solution. We state and prove the theorem only for 2 colours, however, 
the statement as well as the proof are valid for any number of colours. 
Define for / in TC^ 

<M t ,f> = <Y t ,f>-<Y ,f>-[ < l -D 2 (pTY s ,/\f>ds. (6.12) 



We also define the operators S tt i = \j ipi(pt)V, A\ = y/ S'(p t )pt,iA and Vf u the 
semigroup associated to the operator ^S(p t )A. 

Theorem 6.7. Fix a positive integer m > 2. Let P be a probability measure on the 
space {(C[0, T], H^_ m ), J-}, where T = Ut>oJ 7 t, an d Tt is the canonical filtration of 
the process < Y t , f >. Assume that for any f E C°°(T) x C°°(T) and ol e R 2 the 
quantities < M t , f > defined in and 

{a- < M t , f >} 2 - Ei=i l2 «i Si W&Mlds (6.13) 

are L l (P) continuous martingales with respect to Tt. Then, for any colour i, for all 
< s <t, f in C°°(T), and subsets A ofR, 



P[|< Y tti ,f > -J < Y u ,A l u V s (t,u)f > du\ e A\Fs 



(y- <Y i;S ,V t s J >) 
ex P s - „ rt„~ ~c „,o : f d Vi ( 6 - 14 ) 



27T/; \\S u ,{P*J\\ldu I 2 1. \\Su,iV t S uf\\ldu 



P almost surely. 



From this result we obtain uniqueness of the limiting measure P. From the above 
theorem we not only have the distribution of 

<Yu,f>- I <Y u ,AivU> 



for any colour i, but we also know that these quantities are independent for different 
i and different test functions / because of ()6.13|) . If we choose the same / and sum 
over the colours we obtain 

-t 



<Y t ,f>- I <Y u ,A u V t s J>, 

J s 
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where A u is the sum over A l u . This fact plus Theorem 15.81 which gives the full 
distribution of the process < Y t , f > is enough to obtain uniqueness of the pair 
{< Y^, f >,< Yf,h >}. A standard Markov argument gives the uniqueness of the 
finite dimensional distributions, which implies the uniqueness of the measure P. 

The proof of the above theorem is similar to the proof of the uniqueness result for 
the colour-blind density fluctuation field as it appears in [KJL or HSJ. 

Proof. First of all notice that the martingale < M t ^, f > in 1)6.12)1 may also be written 

as 

< M a ,i,f> + < ZlJ > - [ < Y u ,A l J > du, 

J s 

where < Z\ s , f > is used to denote 



< Y tii , f>-< Y s>i , f > 



J < Y Ut i,^S(p u )Af > du. 



Ito's formula together with the fact that (J6.12)) and (j6.13|) are both martingales 
imply that for each fixed s > and / G C°°(T), X tiS (f) defined by 

X tiS (/) = exp|i^<4 s ,/>-^ <Y u ,Aif>di?j+^^ \\S u>i f\\ldu 

is a martingale. Next take two fixed times t\ < ^ < T and let s n j denote a partition 
of the time interval defined by s n j = t± + 1 (t 2 — ti). Consider the quantity defined as 

n" = o-^ s „ iJ+ i,s nj j (V s (T, s n j)f). 
By continuity of < Y*,V s (T,t)f > in {s,t} we may let n —>■ oo to obtain that the 
above quantity converges a.s. and in ^(P) to z * 2 ^ , where Z t (f) is equal to 

exp{i<Y t „Vy>- jf <Y u ,A l u V^J>du+ l -j^ \\S uA V T , u f\\lduY 

Because the convergence above takes place also in ^(P) the martingale properties 
of X t s (V^ Ts ^f) are passed onto Z t (f). The fact that Z t (f) is a martingale proves the 
theorem. □ 

6.4. Tightness. Let P N denote the probability measure on -D([0, T], 7{L m ) induced 
by the fluctuation field Y t N . 

For a function F = {F 1 ,^ 2 } in -D([0, T], 7{L m ), define the (uniform) modulus of 
continuity, Us(F), for a fixed 5 > 0: 

u s (F) = sup {\\F? - F^W-m + \\F? - F s 2 ||- m } . 

\s~t\<S,0<s,t<T 

As previously we shall simplify the supremum in the above notation to sup si . By 
Prohorov's theorem and Arzela-Ascoli it follows that to show the sequence Pjv is tight 
we need to show that it satisfies two conditions: 
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(Tl) Hindoo lim sup N ^ OQ P N [sup < t < T | \Y t \ | _ TO > B] = 
(T2) lim^ limsup^^ P N [u s {Y) > e] = Ve > 0. 

Proposition 6.8. The sequence of measures P^ is tight in D([0, T], Ti.'L^). Moreover, 
all limit points are concentrated on continuous paths. 

Using ()6.2j) we have that 

< Y t N 1 J> = <Y N 1 J> + f < Y s N 1 ,A N f > ds+ < M t N 1 J> 1 

Jo 

where Mf{ is a martingale. By assumption (F2) there is nothing to prove for the 
initial field. Also, M£[ has quadratic variation given in (|6.11|) with a\ — 1 and 
«2 = 0. This formula is very similar to the one in the single colour case, and hence, 
the proof will follow as before from the work we have already done in the two-colour 
Boltzmann-Gibbs principle. We may also handle each of Y/f separately. 
To show that Y t N satisfies condition (T2) it is enough to show 

(T2) A lim (5 ^olimsup iV ^ 00 Pjv[sup Sit sup|| / | U < 1 {< Atv / ', J*Y^' % du >} > f] 
(T2) B lim^ohmsupjv^Pjvtsup^HM^ ~ M s,i\\-m > f], 
for each of % = 1, 2. 

We begin with (T2)b, and fix i — 1. Clearly, the other case is the same. Because of 
this special "decoupling" that occurs, the proof is now almost identical to the colour- 
blind case, and we omit many of the details. The next three lemmas follow from the 
quadratic variation result in Proposition 16. 61 as before; see the proof of Lemmas 15. Ill 
EH21andlo~T51 

Lemma 6.9 (Bounding the tails 1.). There exists a finite constant C = C((p) such 
that for every eigenfunction e z ofTL m , 

limsupPjvf sup | < M^[,e z > | 2 ] < CT{1+ < A N e z ,A N e z >} 

N^oo 0<t<T 

Lemma 6.10 (Bounding the tails 2.). For m > A, 




Lemma 6.11. For any f in C 3 (T) 
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These three lemma together prove that condition (T2) b is satisfied. We next turn 
to condition (T2)a' yet again we split this field up in the sum of several terms: 



= 






*2? = 


s/N{c Kjl (r,?(x)) - ^{mf{x))}l[mf{x) < 


R] 


= 


VNd llfl (p t (x)){m« t (x) - Pu {x)}l[mf{x) 


< R] 


*ff = 


^Nd 2 i Pl (p t (x)){m^(x) - p 2>t (x)}l[mf(x) 


< R] 


= 


y/N [^i(mf(x)) - fi( Pt (x)) 





(6.15) 



-J2diipi(p t (x)){m^ t (x) - pi, t (x)} 

i=l 

S 1 = -VNtp 1 {p t (x))I[mK(x)>T] 
>/Nc Ktl (x)I[m? (x) > r}. 

We need to prove that 



l[mf (x) < R] 



lim lim sup 



sup sup 

s,t II/IU<1 



> 



28 



(6.16) 



for i = 1, . . . 7. As before, the quantities i = 2,3,4 are more difficult. 

For i = 1, 5, 6, 7 we have that 1)6.16)1 holds because we have already shown in the 
previous section 



lim lim E 

1^0 V — >oo 



N 



sup 



r <<'\a^> 

»<1 ./0 



0. 



That is, the result follows by summation by parts for % = 1, and by Theorem 16.41 and 
I ^^0*0 1 < C J2i=i 2^ m f ( x ) — P l ( x )) 2 f° r the remaining cases. 

It remains to show ()6.16jl for i = 2,3,4. By Propositions 15 . 1 51 and IBTHl this follows 
if we have the following bound for i = 2, 3, 4 



E 



EQl 



exp{aN^ 2 (t - s)- 1/2 < A N f, ^ *»Ju >} 



< e 



Ca< 



+1} 



for a positive constant C and all large enough N. As before, to obtain convergence in 
m > 4 we need k = 3. Indeed this follows from Lemma l6~5l and the Sobolev inequality 
for i = 2. It is straightforward to check that the exact same bounds apply for i = 3,4. 
The constant C depends on I and is valid for all a > a , that is, C = C(l,a ). This 
completes the proof of (T2)a- 
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To complete the proof of tightness we hence need the bound (Tl). Because of 
Lemma [6.91 and lb' . 1 Ul we have that for i — 1,2 



We need only 



lim lim sup P/v 



lim lim sup P/v 



sup ||M^||_ m > M 
0<t<T 



0. 



Y N \\ > M 

1 t,i \ \-m 







0. 



again for i = 1,2. However, this follows from Proposition 15 . 1 7l and Lemma f6. 51 This 
completes the proof of tightness in the colour case. 
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